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(U)  The  SEATIDE  Analysis  Process  is  a semi-automated 
procedure  for  the  generation  of  time-phased,  high  value  cruise  missile 
weapon  systems  concepts,  together  with  the  supporting  technology  and 
intelligence  indicators  which  would  reflect  that  these  technological  goals 
are  being  achieved.  The  SEATIDE  process  can  also  be  used  to  evaluate 
the  effectiveness  of  fixed  force  levels,  existing  forces  in  SAL  environ- 
ments, or  Naval  defenses. 

^ (U)  The  Defense  Intelligence  Agency,  through  its  Directorate 

I of  Estimates,  and  The  Advanced  Research  Projects  Agency  (ARPA)  have 

sponsored  the  development  of  this  computer  based  analysis  at  the  weapon 
system  and  Naval  force  structure  level.  A previous  process,  RIPTIDE, 
was  developed  for  DIA  for  use  in  analysis  of  strategic  missile  systems. 

(U)  Generic  to  the  SEATIDE  Analysis  Process  are  three 
major  computer  models;  The  Naval  Engagement  Model  (NEM),  Cruise 
Missile  Concept  Generation  and  Screening  Model  (CM-CGSM)  and  Relative 
Worth  Model  (RWM).  The  NEM  evaluates  force  effectiveness,  tactics,  and 
task  force  configurations;  the  CM-CGSM  enables  definition  and  selection 
of  candidate,  advanced  cruise  missile  system  concepts;  and  the  RWM  per- 
mits assessment  of  worth  in  accordance  with  a variety  of  objective  and 
subjective  criteria.  Each  of  these  models  has  been  checked  out  by  DIA. 

(U)  In  addition  to  exercising  the  computer  models,  there  are 
several  other  analytical  and  engineering  tasks  to  be  performed,  e.g.  , the 
identification  of  areas  of  current  interest  and  the  associated  criteria  and 
potential  concepts,  the  creation  of  a foreign  technology  data  bank  in  a 
format  needed  by  the  computer  models,  the  engineering  of  concepts  to 
the  required  detail,  and  the  use  of  a verification  analysis  loop. 
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AERODYNAMICS  MODEL 


1.0  INTRODUCTION 

The  Aerodynamics  Model  is  designed  to  provide  the  aerodynamic 
parameters  required  by  the  CM-CGSM  for  the  general  solution  to  the  cruise 
missile  synthesis  problem.  The  coefficients  required  for  performance 
analysis  (C  , C ) are  the  primary  output;  however,  certain  elements  of 
the  pitch  and  roll  characteristics  of  the  configuration  are  also  computed  in 
the  module.  The  Missile  Synthesis  ajid  Performance  Routine  (Ref.  1)  was 
used  extensively  as  a model  for  the  SEATIDE  Aerodynamics  Model. 

2.0  OBJECTIVE 

The  objective  of  the  program  is  to  provide  a description  of  the 
basic  force  coefficients  and  pitch  plane  aerodynamic  derivatives  as  a func- 
tion of  vehicle  geometry,  flight  conditions,  and  center  of  gravity  location. 

Force  coefficients  (C.  , C ) required  for  mission  performance  analyses 

JLi  iJ 

are  provided  to  the  Vehicle  Performance  Model. 

3.  0 MODEL  DESCRIPTION 

The  model  computes  aerodynamic  coefficients  of  each  element  of 
the  configuration  and  combines  these  elements  along  with  appropriate  inter- 
ference factors  into  a description  of  the  complete  configuration.  Component 
drag  coefficient  is  computed  as  a function  of  both  Mach  number  and  altitude, 
because  of  the  dependency  of  skin  friction  drag  or  altitude  while  lift  curve  slope 
is  computed  as  a function  of  Mach  number  only.  The  general  procedure  for 
this  is  diagrammed  in  Figure  1.  The  analysis  is  limited  to  the  linear 
angle -of -attack  region,  Mach  numbers  from  0 to  5.  0 and  altitudes  below 
100,  000  feet. 

3.  1 Configuration  Options 

The  configuration  elements  available  for  use  in  the  SEATIDE 
CM-CGSM  are  shown  in  Figures  2,  3,  4,  and  5.  Body  options  are  illus- 
trated in  Figure  2 where  the  various  combinations  of  nose  shape  and  boat- 
tail  are  shown.  The  body  must  have  a cylindrical  center  section  with  one 
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TRAPEZOIDAL  DELTA 


STRAIGHT  TRAILING  EDGE 


TRAILING  EDGE  SWEPT  AFT 


FIGURE  3 LIFTING  SURFACE  PLANFORM  OPTIONS 
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tail  arrangements 


PLANAR  TRIFORM 

CRUCIFORM  + CRUCIFORM  X 


SECTION  PROFILES 


DOUBLE  WEDGE  BI  CONVEX  MODIFIED  DOUBLE 

WEDGE 


FIGURE  4 LIFTING  SURFACE  ARRANGEMENTS 
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BELLY  MOUNTED  INLETS 


LAUNCH  LUG 


CAPTURE 

AREA 


SIDE  MOUNTED  INLETS 


COWL  LIP 


PROJECTED  AREA 


inlet  capture 

AREA 


BOUNDARY 

LAYER 

DIVERTER 

WEDGE 


AFT  FAIRING 


FIGURE  5 CONFIGURATION  OPTIONS  - INLET  SYSTEMS 
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TABLE  I 


AERODYNAMIC  OPTIONS 


NAME 

DESCRIPTION 

OPTION 

NO 

ITN 

Type  of  nose 

Tangent  Ogive 

1 

Von  Karman 

2 

Cone 

3 

Spherical 

4 

Blunted  Cone 

5 

Blunted  Ogive 

6 

IBTL 

Type  of  boattail 

None 

7 

Conical 

8 

ILUG 

Launch  lugs 

No 

9 

Yes 

10 

NW 

Fixed  wings 

No 

11 

Yes 

12 

ICNTRL 

Control  Selection 

Tail 

13 

lART 

Aerodynamic  Surface 

Planar 

14 

Arrangement  (Aft) 

Cruciform  + 

15 

Triform 

16 

Cruciform  X 

17 

lARW 

Aerodynamic  Surface 

Planar 

18 

Arrangement  (Fwd) 

Cruciform  X 

19 

ITSE  CT 

Aerodynamic  Surface 

Double  Wedge 

20 

Section  (aft) 

Bi-Convex 

21 

Modified  Double 
Wedge 

22 

iwsect 

Aerodynamic  Surface 

Double  Wedge 

23 

Section  (Fwd) 

Bi- Convex 

24 

Modified  Double 
Wedge 

25 
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TABLE  I (Cont'd.  ) 


NAME 

DESCRIPTION 

OPTION 

NO. 

ISURFW 

Type  of  wing  planform 

Trapezoidal 

26 

Delta 

27 

ISURFT 

Type  of  tail  planform 

Trape  zoidal 

28 

Delta 

29 

IPLANW 

Wing  planform  definition 

Input;  SEW,  ARW, 

30 

input  options 

TRW,  = 0.  0 

’ TE 

Input:  SEW,  ARW, 

\ , \ =0.0 
•'LE’  TE 

31 

Input:  SEW,  ARW, 

32 

^DES’  TE  " ° 
Input;  SEW,  ARW, 
TRW,  A,t,^ 

33 

IPLANT 

Tail  planform  defini- 

Input; SET,  ART, 

34 

tion  input  options 

TRT,  =0.0 

Input:  SET,  ART, 

\ . \ =0.0 
LE’  ‘'TE 

35 

Input;  SET,  ART, 

36 

M^^^,  A ,^^  = 0.  0 
DES  TE 

Input;  SET,  ART, 

37 

TRW, 

TE 


r 


of  the  nose  shapes  shown,  and  may  have  a conical  boattail.  Nose  fineness 
ratio,  bluntness  ratio,  and  boattail  fineness  ratio  are  input  variables.  The 
amount  of  boattailing  desired  may  be  selected  by  an  input  factor  which  de-  | 

fines  the  amount  of  the  base  annulus  to  be  removed  by  boattailing.  | 

Lifting  surface  options  are  illustrated  in  Figures  3 and  4.  Lift- 
ing surfaces  may  include  two  sets  of  surfaces  (wing-tail)  or  a body-tail 
configuration.  Arrangement  of  the  surfaces  may  be  planar  (2  panels), 
triform  (3  pzuiels),  cruciform  +,  or  cruciform  X as  shown  in  Figure  4. 

Wing  and  tail  panels  are  assumed  to  be  in  line  in  this  model. 

Trapezoidal  or  triangular  planforms  may  be  specified  by  input.  Wing  plan- 
form  characteristics  are  specified  in  all  cases  by  exposed  area  and  aspect 
ratio.  Taper  ratio,  leading  edge  sweep,  Mach  number  for  a sonic  leading 
edge,  or  trailing  edge  sweep  angle  may  be  used  in  various  combinations  to 
select  a surface  planform.  Table  I shows  the  input  options  available  for 
surface  planform  definition.  Figure  6 illustrates  the  lifting  surface  plan- 
form  parameters. 

Option  30  requires  input  values  for  exposed  wing  area  (SEW), 
exposed  aspect  ratio  (ARW),  and  taper  ratio  of  the  exposed  panel  (TRW); 
the  trailing  edge  is  unswept.  Exposed  span  (BW)  of  the  wing  is  defined 
from  the  relationship  between  aspect  ratio,  span,  and  area; 


ARW 


BW^ 

SEW 


BW  = ^SEW  =:=  ARW 


Wing  root  chord  (ROW)  is  defined  by  the  relation 
ROW  = 2 X SEW/BW/(1  + TRW) 

Tip  chord  is  defined  by  the  taper  ratio  and  root  chord 


TCW  = ROW  (TRW) 


I 
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THEO 


FIGURE  6 


CONFIGURATION  PARAMETERS 


Leading  edge  sweep  angle  is  defined  by  the  relation 


JV  = TAN-‘  ^ 

•''■LE  ARW  (1.  + TRW) 

Option  31  requires  input  values  for  leading  edge  sweep  angle 

rather  than  taper  ratio;  the  trailing  edge  is  unswept.  Taper  ratio  is 

computed  from  the  following 


ARW  TAN-^LE 


TRW  = 


) 


1.  + ARW  TAN-^LE 


If  the  input  leading  edge  sweep  combined  with  the  input  aspect 
ratio  is  not  consistent  with  a trapezoidal  wing,  a delta  wing  with  the  input 
leading  edge  sweep  is  defined  from  the  input  area  and  a computed  value  of 
aspect  ratio. 

ARW  = 4.0/TANA  _ 

LE 

BW  = \SW  ARW 

Wing  root  chord  for  the  delta  wing  is  based  on  the  relationship 
ROW  = 2.  0 * SEW/BW 
Tip  chord  is  zero  for  the  delta  wing. 

Option  32  is  similar  to  Option  31  with  the  leading  edge  sweep  angle 
being  set  to  provide  a sonic  leading  edge  at  a specified  Mach  number.  Lead- 
ing edge  sweep  angle  for  a sonic  leading  edge  is  defined  by  the  relation 

A le  = tan"^  Vm^  - 1 


The  remainder  of  the  wing  parameters  are  computed  as  in  Option  31.  Trail- 
ing edge  sweep  angle  must  be  input  in  Option  33  along  with  aspect  ratio, 
exposed  area,  and  taper  ratio.  Leading  edge  sweep  angle  is  computed  from 
the  relation 


- 1 B W t 

Ale'"’^^^  (2(RCW+^  TAnA.j.^  - TCW)/BW) 

Midchord  sweep  angle  is  computed  from  the  following 
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A =tan’^  RCW-TCW  a 

Quarter  chord  sweep  angle  is  based  on  the  relation 


JV  , --  TAN-‘  ■■M.R^.-.TCWg  ,71 

C/4  BW  TE 

Mean  geometric  chord  length  is  based  on  the  relation 

2 

MGC  = 3 ( ^ 

Spanwise  location  of  the  mean  geometric  chord  is 


_ ^ /I.  + 2 TRW  \ 

^MGC  ~ 6 \1-  + TRW  ) 

Longitudinal  location  of  the  leading  edge  of  the  mean  geometric  chord 
with  respect  to  the  surface  leading  edge  intersection  with  the  body  is 

^MGL  ~ ^MGC 

The  location  of  the  intersection  of  the  wing  leading  edge  may  be 
input  or  computed  to  satisfy  a desired  stability  margin  at  either  a subsonic 
condition  (M  = 0.  8)  or  a supersonic  condition  (M  = 2.0)  for  example. 

Section  profile  may  be  double  wedge,  bi-convex,  or  modified 
double  wedge  as  shown  in  Figure  4.  Surface  thickness  ratio  is  an  input 
variable. 

Two  inlet  configurations  are  provided  for  airbreathing  configu- 
rations as  shown  in  Figure  5 . A single  two-dimensional  belly  mounted 
inlet  or  dual  side  mounted  inlets  may  be  selected.  Launch  lugs  may  be 
specified  as  required. 

3.  2 Input -Output  Options 

Input  requirements  are  related  to  the  configuration  options 
selected.  The  basic  inputs  are  body  diameter,  missile  length  or  weight, 
wing  area,  wing  aspect  ratio,  and  tail  area.  Modifying  inputs  for  the  body 
are  nose  shape,  fineness  ratio,  and  bluntness  ratio.  Wing  and  tail  modifi- 
cations include  sweep  angle,  taper  ratio,  section  profile,  and  thickness 
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Table 


ratio.  Inlet  type  must  be  specified  for  the  airbreathing  configuration. 

I lists  the  configuration  option  switches  for  the  Aerodynamics  Model.  De- 
tailed input  instructions  are  discussed  in  Volume  IIIA. 

The  basic  output  of  the  aerodynamics  model  is  zero-lift  drag 
coefficient  and  lift  curve  slope  of  a specified  configuration  as  a function  of 
Mach  number  and  altitude.  Table  II  lists  the  aerodynamic  coefficients 
computed  in  aerodynamics  model  and  included  in  the  lift/drag  coefficients. 
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TABLE  II 


AERODYNAMIC  COEFFICIENTS 


Symbol  Description 

CdO^t^^  Zero-lift  drag  coefficient  (power-off) 

OFF 


cdo 


ON 


Zero-lift  drag  coefficient  (power-on) 


m , 


CmQ 

^rn 


Cl 


P 


Ch 


d( 


Slope  of  lift  coefficient  curve  with  respect  to 
angle  of  attack  d C 

j-i 

d cX 

Slope  of  lift  coefficient  curve  with  respect  to 
surface  deflection  d C 

Lj 

Slope  of  pitching  moment  curve  with  respect  to 
angle  of  attack  d C 

m 

d o( 

Slope  of  pitching  moment  curve  with  respect  to  surface 
deflection  d C 

m 


Slope  of  pitch  damping  moment  with  respect  to 
pitch  rate  d 

Slope  of  pitch  damping  moment  with  respect  to  rate  of 
change  of  angle  of  attack  d 


Slope  of  rolling  moment  coefficient  with  respect  to 
surface  deflection  d 


Slope  of  roll  damping  moment  with  respect  to  roll 
rate  d 


Slope  of  tail  hinge  moment  curve  with  respect  to  angle 
of  attack  d (HM) 
d o<L 


A-18 


TABLE  II  (Cont'd.) 


Symbol 


CNr 


c 


Y 


Description 

Slope  of  tail  hinge  moment  curve  with  respect  to  tail 
surface  deflection  d (HM) 

d 5 

Slope  of  yaw  damping  moment  with  respect  to  yaw 
rate  d 

dR 

Slope  of  yawing  moment  with  respect  to  sideslip 
angle  d 

d B 

Slope  of  side  force  coefficient  with  respect  to  sideslip 
angle  d 

Ta 
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3.  3 Analysis 

Complete  configuration  aerodynamic  coefficients  are  deter- 
mined by  a build-up  procedure  with  appropriate  interference  factors.  The 
general  procedure  for  this  is  diagrammed  in  Figure  1 using  the  component 
options  described  in  Table  I. 

The  procedures  used  for  computing  the  component  coefficients 
and  interference  factors  are  outlined  in  the  following  paragraphs,  along  with 
the  source  references. 

4.  0 Body  Aerodynamics  Characteristics 

The  vehicle  body  consists  of  a right  circular  cylinder  with  a 
forebody  and  possibly  a conical  afterbody.  Forebody  shape  and  fineness 
ratio  are  input  variables  while  afterbody  fineness  ratio  and  diameter  ratio 
are  variable. 

4.  1 Drag 

Body  drag  consists  of  skin  friction  drag,  nose  pressure  drag, 
afterbody  pressure  drag,  base  drag,  and  drag  due  to  lift. 

4.1.1  Zero-lift  Drag 

Subsonic 


L CFRBy- 


+ 


ll 


SH£.~r 

SgSP 


The  flat  plate  incompressible  skin  friction  coefficient  (Cp  J is  based  on  the 

modified  Prandtl-Schlicting  relationship  for  turbulent  boundary  layers. 

_ a.A8Z 

Body  Reynold's  number  is  computed  from  1962  U.  S.  Standard 
Atmosphere  data.  Reynold's  number  per  foot  per  Mach  number  data  vs. 
altitude  is  presented  in  Table  III.  The  incompressible  skin  friction  coef- 
ficient is  corrected  for  compressibility  effects  by  the  data  of  Table  IV. 

Base  drag  at  all  Mach  numbers  is  computed  from  the  base 
pressure  coefficients  presented  in  Table  V.  Base  drag  coefficients  are 
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TABLE  III 

REYNC'LD'S  NUMBER  PER  E'i  . PER  MACH  NUMBER 

VS  ALTITUDE 


data 

F■KCT^t/  r. 

, 7.0P>n, 

c; 

• 

, ^.171, 

1 r-  . 

IS. 

, '-.''11. 

★ 

?0  . , 

?S. . 

■^.370, 

3^'., 

r 

m 

? . '1  r 1 , 

★ 

ilO.  , 

1 .“51  u. 

! . 1 

.'’71, 

* 

^0. , 

.7?<5, 

, , 

.^7k. 

. a ^ , 

7-^.  . 

.IS'- , 

* 

PO., 

.?Pn, 

'3'^. , 

.171, 

'5'=.  , 

. 1 7A, 

* 

1 '^'1. , 

. 1 ''7  / 

TABLE  IV 

COMPRESSIBILITY  CORRECTION  FACTOR  VS  MACH 
NUMBER  AND  ALTITUDE 


niTA  FPCT^l3/ 

O.n, 

c 

• 

c 

1 .r*. 

o.q7-5. 

f . 7 . 

1 

0 .ouS, 

1 .fl. 

n.pp-. 

1 .S, 

0 . , 

7.", 

2 

3.^, 

n .S7^j , 

a. 

'•■  . u 'J  S . 

3 

r .3‘'3, 

'-.'1. 

f’.pp. 

«r.o. 

r..n. 

1 

a 

n.s. 

o.o7«. 

f’  .7t  , 

1 .r'. 

o.q37. 

5 

1 .s, 

O.BSS, 

?.s. 

0 . f-  . 

3.r., 

''.‘51 . 

h 

<J  , 0 . 

'i.aqi  , 

S.r , 

r. . a , 

'-.r . 

^ . 3 a 7 / 

TABLE  V 

BASE  PRESSURE  COEFFICIENTS  VS  MACH  NUMBER 
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corrected  for  boattailing  by  the  base  diameter  ratio. 

Power-on  base  drag  coefficient  is  computed  from  the  relation: 


Power-off  base  drag  coefficient  is  based  on  the  relation: 

Transonic,  M=1 

Nose  pressure  drag  at  M = 1 is  based  on  curve  fits  of  the  data  of 

Reference  3.  Tables  II  and  III  present  ogival  nose  and  conical  nose  pressure 

drag  coefficients  at  M = 1 as  a function  of  nose  apex  half  angle.  Boattail 

pressure  drag  coefficient  at  M = 1 is  based  on  data  of  Reference  4.  Table  III 

presents  boattail  drag  parameter  _ /£.  Osr  as  a function  of 

J & 

boattail  base  area  ratio.  The  effect  of  nose  bluntness  on  nose  drag  coef- 
ficient at  M=1  is  based  on  data  of  Reference  5.  Base  drag  and  friction 
drag  coefficients  are  computed  in  the  same  manner  as  the  subsonic  condi- 
tion. Launch  lug  drag  coefficient  at  M=1  is  based  on  empirical  data  of 
Reference  6.  Zero-lift  drag  at  M=1  is  given  by  the  relation: 


Supersonic,  M = 1.  25  - 5 

Nose  pressure  drag  coefficients  for  tangent  ogive,  von  Karman, 
and  conical  nose  shapes  are  based  on  a curve  fit  of  the  data  of  Reference  7. 
Tables  IX  through  XI  present  nose  pressure  drag  coefficient  times  the  square 
of  Mach  number  as  a function  of  Mach  number  divided  by  nose-fineness  ratio. 
Hemispherical  nose  pressure  drag  coefficient  from  Reference  8 is  presented 
in  Table  XII  as  a function  of  Mach  number.  Nose  pressure  drag  coefficient 
of  hemispherically  blunted  conical  noses  is  based  on  data  of  Reference  1. 
Tabular  data  of  nose  pressure  drag  for  blunted  cones  is  presented  in  Table 
XIII  as  a function  of  Mach  number,  theoretical  nose  fineness  ratio,  and  nose 
bluntness  ratio.  Pressure  drag  coefficients  for  hemispherically  bl’-nted 
ogive  noses  is  presented  in  Table  XIV,  as  a function  of  Mach  number,  theo- 
retical nose  fineness  ratio,  and  nose  bluntness  ratio. 
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2 
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1 

2 
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O.P70 , 

70.000, 
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00.000, 
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TABLE  VIII 
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r.-PK  T, 

on? 
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0 .2P0  , 
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0.110, 
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*-AC1' 

S(ftlJA9 

CPF 

D’^L 

cnwt-'2 

DF'l 

CP1'”2 

PAJA  COmVK/  O.OOO, 

0.0(10, 
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0 , 8 9 ft  , 
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0 . 1 P 0 , 

3 
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u 

'1.317, 
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5 
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O.PFC, 
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6 
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7 
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P 

1 
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<J 
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X 
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r 
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X 
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0.170, 
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0 
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X 
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0.070,  0,045. 
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0 
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X 
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ft 
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X 
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0.152. 
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X 
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1.151  , 
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X 
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ft 
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O.PO 
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1,00 

3 
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0.202,  0.2‘<<. 

0.711, 
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ft 
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a 
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0.42P. 

ft 
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5 
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0.245, 
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ft 
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b 

0 . 0 ® . 
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0.232. 
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ft 
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7 
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r 
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P 

0 . 0?F  , 
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ft 
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9 
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O.P‘JS, 
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F 
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G 
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P.2P6, 

0 .308  . 

0 . il  7 8 . 

0.7CO , 

R 

O.OPS, 

P.  1 OP, 

0.183. 

0.281 , 

0 , ii  7 3 , 

0.753, 

c 

n.CRS. 

0.068, 

0.  1-^8, 

0.269, 

0.868, 

''.75iJ  , 

D 

n.03S, 

P.PG?, 

0.  12b, 

0.263, 

0.869, 

0.75-3, 

F 

P.P25, 

P . Pa? , 

0.1?a, 

0.260, 

0.a79, 

0 . 75iJ  . 

F 

P.018, 

P.P30, 

0.115. 

P.?87, 

P . a88  , 

P.75il, 

G 

O.OIS, 

P.P36, 

0.115. 

0.25P, 

p . a 6 g , 

0.75U, 

PWS^.PO, 

FPs?.p  T8RI 

8.0 

FCAPsP.P 

0.2  0 

P.aP 

0 ,hP 

P.8  0 

1 .OP 

X 

0.173, 

P.178, 

0.?‘^8, 

P.36P, 

6.S8a, 

0.8 30, 

1 

0.088, 

0.111. 

0.185, 

0.33}  , 

p.aag. 

6 . P30, 

2 

0.0'^3, 

P . PP  P , 

0.168, 

0.323, 

P . 5 U 1 , 

P . a t 0 , 

3 

0.036, 

P . P68 , 

0.156, 

P . 31 3, 

P.638, 

o.5a6. 

<i 

0.026, 

P , P88 , 

0.180, 

P.31'3, 

P.53P. 

6 . a <4 1 . 

P.01  «>, 

P .P60  . 

P . 1 ab  , 

0.3(,R, 

.5  38  , 

p,8a?. 

h 

P.P1  a. 

P.  Pii", 

0.102, 

0.366, 

0.537, 

6 . at.  1 / 

P8sa , PP  , 

FPsP.P 

8.0 

FCAPsP , P 

0.20 

P .UP 

0/  . 6 P 

6,8'J 

1 . 0 P 

pata 

rPp713/  0.163, 

0.173, 

0 .?ao , 

0.38,6  , 

0.800, 

0.868, 

8 

''.081  , 

P . 1 1 1 . 

0.20'-, 

P.355, 

0.577, 

0.567, 

0 

''.OFP  , 

''.  08  2 . 

r . 1 80  , 

0.300 , 

0.588  , 

0.86  8 , 

A 

P.P32. 

0 . 0 7 0 , 

P . 1 7^  , 

P.338, 

0.585. 

0.860, 

H 

0.P22, 

P . 0 r-  ^ , 

0.182, 

0.333, 

0.581  , 

0.860, 

c 

P . 0 1 6 , 

P . 0 8 a , 

0.18,-,, 

P.33P, 

0.560, 

0.869, 

D 

0.011, 

P.P88, 

0.186, 

0.3^5, 

6.860, 

0.570/ 

y?  5 ' • 

■■  ^ i r 
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Boattail  pressure  drag  is  based  on  data  of  Referent  e 1.  Table 
XVpresents  boattail  drag  parameter  as  a function  of  boat- 

tail  base  area  ratio  ( boattail  parameter 

Base  drag  is  computed  as  described  in  the  subsonic  region. 
Body  friction  drag  is  based  on  the  relation 

<^F<.  Sit.ep 

Launch  lug  drag  coefficient  is  based  on  empirical  data  of 
Reference  6.  Table  XVI  presents  launch  lug  equivalent  flat  plate  area  as  a 
function  of  Mach  number.  Launch  lug  drag  coefficient  is  based  on  the 
relation  ^ 

C = 

Sgs  r 

Zero-lift  drag  coefficient  at  supersonic  Mach  number  is  based 
on  the  relation 


■£>c 


= c 


-Oar 


C 


oa 


4.  1.  2 


relation 


4.  1.  3 


Induced  Drag 

Induced  drag  foeff i<'ient  at  all  Mach  numbers  is  based  on  the 


Total  Drag  Coefficient 
Total  drag  coefficient  is  based  on  the  relation 
Co  ^ +Coi. 

4.  2 Lift  Curve  Slope 

Body  lift  consists  of  the  nose  contribution  and  the  boattail 
contribution. 

Subsonic,  M 9 

Nose  and  cylindrical  section  lift  curve  shape  is  based  on  poten- 
tial theory  and  uses  Munk's  apparent  mass  factor  from  Reference  1. 

Munk's  efficiency  factor  is  presented  in  Table  XVII  as  a function  of  body  fine- 
ness ratio.  The  boattail  contribution  is  equivalent  to  the  nose  contribution 
corrected  for  boattail  volume  and  boattail  base  area  ratio.  The  boattail 
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TABLE  XV 


CnViICAI.  ”f'4TTAII  rt,tsshl.>  rp.‘.r;  OAWflMt-Tf;': 

''TPstf^A'^p  PTA./  •'hi , 

CPAsCPPBTl  /nw*p 

W1  OHsMTiA  TTA  Tl  PAPA^FJPP  f ? , n / J F T A ) * f ( / P 1 


PI 

DPsI  . 0 

data 

CPABTI / 

1 . p 0 n , 

r-n^ 

PDA 

r.',p 

r- A 

rr? 

r''A 

1 

• 

n.oon. 

0 . ap  0 , 

0.025, 

0.570, 

0 . 0 5 n , 

O.50O, 

2 

n.f)7S, 

0 .5?o  , 

0 . 1 0 0 , 

0 . 5 0 0 , 

0.  150, 

0.5fiO, 

3 

P,2on, 

0.51 0, 

o.2':o. 

0 . aHo , 

0.300, 

0 , <i  'J  5 , 

a 

^.3F0, 

0. 3P0, 

0 . a 0 0 , 

0.7?0, 

o.a5A, 

'1  . P50  , 

5 

0.500, 

0.?1 0, 

O.550, 

0 . 1 50 , 

O.ooo, 

0.1 30, 

6 

0,R50, 

0 . 1.  1 0 , 

0.700, 

0 . 090  , 

0 . F 0 0 . 

0.050, 

7 

0.B5O, 

0 . oao , 

0.'>5o, 

0.000, 

1.000, 

,0.000  , 

DPs?.  0 

R 

2.000. 

DP? 

CPA 

P['P 

ri'A 

I'i'? 

r PA 

q 

0.000, 

1 .7ao  , 

n.0?5. 

1 .5ao, 

0 . 050  , 

' . 5 a 0 , 

A 

0.075, 

1 .U/JO  , 

0 .10  0, 

1 .^ao, 

0,150. 

1 .200, 

B 

0 .200  , 

1.070, 

0.250 , 

0 . <0  14  0 , 

0.500, 

^ . F 1 . 

r 

0.35O, 

0.700, 

0 . a 0 r , 

0.590 . 

0 .'-15  0, 

0.500, 

0 

0.500, 

o.a?o. 

0.550, 

0.750  , 

. F 0 0 , 

0.290, 

F 

o.Fsn, 

0.210, 

0.7O0. 

0.170, 

r . p 0 0 , 

O.OfiO, 

F 

0.900, 

0.020, 

0.950, 

0 . 0 1 0 , 

1 . noo , 

0.000/ 

PI 

PP  = 3,0 

DATA 

CPABT?/ 

3.000, 

pp? 

or  A 

r.  r 

COA 

r P? 

roA 

1 

o.ooo. 

2 .URO  , 

f .o?5. 

2.P50, 

0.050, 

2 . opo , 

2 

0.075, 

1.590, 

0 . 1 or  , 

1.750, 

0 . 1 5 0 , 

1 , 5 '1  0 , 

3 

0.200, 

1 . 3ao , 

0.250, 

1.170, 

0.300, 

1.010', 

a 

0.350, 

0.P70, 

0 . a 0 0 , 

o.7?r , 

% 

c 

if 

• 

c 

0.509, 

5 

0.50O, 

0 .50  0 , 

0.550, 

0 . 'M  0 , 

0 . 5 0 0 . 

0.320;, 

A> 

0 . RPO  , 

0.250, 

0 . 7 0 r , 

0 . ? fi  0 , 

0 . fi  r-  0 , 

0 . 0 9 r. , 

7 

''.900  , 

o.o?o. 

0 .950., 

0.O1O, 

1 . OOO , 

0 . 0 r.  0 , 

1 

0 5 = ii  . 0 

fl 

a.ooo  , 

DP? 

r f'  A 

rtv 

fO'A 

r r p 

r oiA 

9 

o.ooo. 

2.0qA, 

o.r?5. 

2.500, 

0.050, 

P.TfiO , 

A 

0.075, 

2.2?0, 

0. 1 or  , 

2.050, 

0.150, 

1 .7P0  , 

e 

0 . 2 0 0 , 

1 .550, 

r . 250  , 

1 . 350  , 

0.300, 

1.170, 

c 

r.^Fo, 

1 . 0 0.  , 

O'  . a 0 r,  , 

O.B50  , 

0.950, 

0.700, 

D 

0.500, 

0 . 7 0 , 

r. . 5 5 n , 

0 . a 5 0 , 

0 . 5 0 r , 

0.750, 

F 

''.550, 

. ?90  , 

f , 7 0 0 , 

0 . ^uO , 

0 . fi  r 0 , 

0 , nq  n , 

F 

0 . 0 0 0 , 

' . 0 ? '•■ , 

0.050, 

0.010, 

1 . 0 0 r , 

0.005/ 

■ • ••'  I 

1 
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SL  np  = «; . ft 

DATA  CDAPT1/  s.r.rr. 


rr? 

rft)A 

or? 

ftOA 

rr? 

r ftA 

ft  , Oftft  , 

7,3?ft, 

ft.ft?5. 

3.ftlft, 

ft . ftpft , 

?,8Sftf 

f).07S. 

? . 'J  ft  ft  * 

ft.lftft. 

?.?70, 

O.l'^O. 

1 .97ft. 

n.2ftc. 

1 .7lft, 

ft,?5ft. 

1 .UBft, 

ft.3ft0. 

1 .?6ft. 

ft.BPO, 

1 .ftBft, 

ft.Uftft, 

ft. 91ft, 

ft.PBft, 

ft  .77ft , 

n.^ftft. 

ft.ft?ft. 

ft.'^'ift, 

ft . p 1 ft , 

ft. ftftft. 

ft.iJl  ft. 

ft.ft'^ft. 

ft. 31ft, 

ft  .7ft0  , 

ft.?Uft, 

ft  .8ft0  , 

ft.  1 1 f'. 

ft  .«»ftft . 

ft  . ft  3 ft  , 

ft.9Sft, 

ft.ft?ft. 

1 .ftftft. 

0,00ft, 

k'l.np 

C 

• 

II 

^).0ftft  , 

f'n? 

ft  fti 

or? 

coi 

Oft? 

CftA 

ft  .ftftft , 

3.7?ft, 

ft . ftpft , 

?,96ft. 

ft.lftft. 

?.97ft. 

0.1‘^ft. 

?.1 1ft, 

ft  .?ftn , 

1 .PUft, 

ft.2‘=0. 

1 .8ftft. 

ft  , 3 ft  ft  , 

1.39ft, 

0.3Pft. 

1 . 1 7ft, 

n.u/'o. 

0 . ORO , 

ft.tipft. 

ft  , Bdft  , 

ft.Bftft, 

ft.7ftft  , 

O.SBft, 

ft  .Plift, 

ft  . ft  ft  0 , 

ft  , M P ft  , 

ft.6?ft. 

ft.3uft. 

ft.7ft0. 

ft.ppft. 

ft.7Pft  , 

ft.lftft, 

ft,  , 8 ft  ft  , 

ft. lift. 

ft.SPft, 

ft  . ft  8 ft  , 

ft.pftft. 

ft,ft3ft. 

ft.'^BO, 

ft.ft?ft. 

1 .ftftft. 

0. ftftft/ 

ci.  msp.ft 


CPAPTil/  P.CftO, 


Oft? 

COA 

00? 

COA 

Onp 

COA 

ft. ftftft. 

u.?5ft. 

ft  .ft?5. 

3.750, 

ft. 0*^0, 

3.?7ft, 

ft.ft7S, 

3.ft?ft, 

ft.lftft, 

?.78ft. 

ft.  I'^ft, 

2.380  , 

ft.Bftft, 

?.ftoft. 

0.?50, 

1.78ft, 

ft. 300  , 

1 .520  , 

ft,3Pft, 

1 .3Cft, 

ft.Pftft, 

1.090, 

ft.9CO, 

0.91ft, 

O.Sftft, 

ft.7feft, 

ft. 550, 

ft. 80ft, 

0.8ftft, 

0.980, 

ft.BPft, 

ft. 37ft, 

0.70ft, 

ft,?7ft. 

0.800, 

0.130, 

ft.9ftO, 

ft  , ftUft  , 

ft  .95ft  , 

ft.ftl  ft. 

1 .nro. 

ft.Oftft, 

PI  OP. 

= 1 ft  .0 

1ft. ftftft. 

Oft? 

ftO« 

00? 

COA 

on? 

roA 

ft. ftftft  , 

5. ftftft  , 

n.o?5. 

9.?7ft, 

0.050, 

3.58ft , 

0.075, 

3.?8ft, 

ft.lftft. 

P.900, 

0.150, 

?.58ft. 

ft.?ftft. 

?.?1 ft. 

ft.?8ft. 

1.91ft, 

0.300, 

1 .‘>30  , 

ft. 35ft, 

1 .39ft, 

0.900, 

1.18ft, 

O.ijPft, 

0.97  ft,- 

ft.Sftft, 

0.p?0, 

ft  .55ft  , 

ft. 880, 

ft  .80ft , 

0.920, 

ft.BPft, 

ft. 91 0, 

P.7ftft, 

0.270, 

ft.Pftft, 

0 . 1 9ft  , 

ft.Pftft, 

ft , ftuft , 

0.950, 

ft.ftl  ft. 

1 .ftftft. 

0 .000/ 

RLDPsI 6. ft 


DATA  Cr)APT5/1h.0ftft, 


or? 

COA 

npp 

C0)A 

r np 

COA 

1 

0. ftftft. 

9.8ft0 , 

ft.0?5. 

P.ftftO, 

0.050, 

9.10ft 

2 

ft. 075, 

3.70ft, 

ft.lftft. 

3.39ft, 

ft.  15ft, 

?.91  0 

3 

ft.POft, 

?.900 , 

0.?9ft, 

?.  1 8ft  , 

0.300, 

1.98ft 

9 

ft. 350, 

1.830, 

ft  . 9 0 C , 

1.380, 

0.990, 

1,19ft 

5 

ft. 500  , 

0.9?0, 

ft. 55ft, 

ft. 750, 

ft  . 8ftft  , 

A , CBO 

8 

0.890. 

0.990, 

0.7ftft, 

ft.  320  , 

o.«ftft. 

0 . 1 90 

7 

0.900, 

0.090, 

0.99ft, 

ft.ftl  *. 

1.00ft, 

0 , rno 
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R I.  n 0 s p n . c 

8 ?0,0OO, 


TP? 

CPA 

PP? 

roA 

PC? 

CO  A 

Q 

a , RRO , 

o.nsn. 

'J  . P 0 , 

0 . 1 0 0 , 

5.7F0, 

A 

0.1SP, 

5.550, 

P.200, 

2.750, 

0 .?‘:0  . 

?.5ft0. 

F 

0,500, 

?. ouo  , 

0.350, 

1 .7ao, 

0 .aoo , 

1 .ii7o. 

C 

O.iJFO, 

1 . ? 0 0 , 

0 . 5 0 0 , 

1.000, 

O.550, 

0 . “00 , 

D 

0 .frOO  , 

0.A50, 

0.550, 

0 . aRo , 

0.700, 

0.550, 

F 

0.7F0, 

O.pto, 

0, , « 0 r , 

0.150, 

0.F50, 

0 , OPO  , 

F 

O.Ron, 

0 . oa" , 

n.<950. 

0.01 0, 

1 , 000  , 

0,000/ 

table  XVI 

LUG  nc«G  t^nPFFTf  IFN,T 

Ft  s LUR  EOUT'^AtFM  PIATF  4RfA-‘?onacf 


RN-a  PACP 

MiimbFP 

BM 

FI. 

0*^ 

FI 

C M 

F| 

DATA  FL0^^7/  1.200, 

18.200, 

1 .250  , 

1 8 . a n 0 , 

1 .30  0 , 

1 0 . a 0 0 , 

1 .iiOO, 

18,100, 

1.600, 

1 6.R00 , 

2. ooo  , 

1 a . Q 0 0 , 

2.PO0  , 

13.200, 

2.800, 

1 1 . 600  , 

3. 20  0 , 

1 n . 1 0 0 , 

3.600, 

9.100, 

U . 0 0 0 , 

8.1 00/ 

table  XVII 

MIIMK  I 5 

FPF  jr 1 F'TV 

F AC  Tf.p 

fph=  ppny 

FT^F^,F5S  RATIf' 

F K s 1 1 M<  • ^ 

FFP  TfTt 

•'TV  FACTOR 

FBB 

Fx 

F RP 

cv 

FPP 

FK 

DATA  FkFPP/  a.oor, , 

. 7P0  , 

5 . 0 0 r , 

o,a?o. 

6 . 0 0 0 , 

O.P70, 

7 . 0 0 0 , 

.poo. 

“ . 0 0 r: , 

o.n^o. 

O.OOO  , 

0.030, 

1 0 . 0 0 n , 

n 

. ®a  0 , 

11.000, 

0.R50  , 

12.000, 

''.opo. 

1 3.0on, 

0 

. '705, 

1 a . 0 n , 

0.070, 

18.000, 

0.075, 

16,000, 

n 

.060  , 

17,000, 

o.oPa, 

1P.O00, 

O.OFF, 

19.000, 

0 

,R«P, 

20.000, 

0 . OP  P / 

factor  is  presented  in  Table  XVIII  as  a function  of  boattail  geometry  and 
angle  of  attack.  Body  lift  curve  slope  at  subsonic  speeds  is  based  on 
the  relation 

- 2 £K  -h  ' 

Transonic,  M=1 

Body  lift  curve  slope  at  Mach  1 is  based  on  slender  body  theory 
adjusted  for  experimental  results.  Body  lift  curve  slope  is  given  by  the 
following  relation 


Supersonic,  M>1.2 

Supersonic  lift  curve  slope  for  nose -cylinder  bodies  is  a func- 
tion of  nose  fineness  ratio  to  cylindrical  section  fineness  ratio. 

Supersonic  lift  curve  slope  for  cone-cylinder  bodies  is  pre- 
sented in  Table  XK.  Ogive-cylinder  data  is  presented  in  Table  XX.  Effects 
of  nose  bluntness  on  blunted  cone  lift  curve  slope  are  based  on  data  of 
Figure  4.  2.  1.  1-24  of  Reference  1.  Tables  XXI  and  XXII  present  bluntness 
effects  for  blunted  conical  nose  shapes  and  blunted  ogival  nose  shapes, 
respectively. 

4.  3 Center  of  Pressure 

Subsonic  and  Transonic 

Center  of  pressure  for  nose  and  boattail  at  subsonic  speeds 
and  transonic  speeds  according  to  potential  flow  theory  is  given  by  the 
following  relation: 


-I- J=§l/fL 


3(0^  - ' ^)7  ^ 
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TABLE  XVIII 

RPATTiTL  lift  FACTrt- 

VIA  RATsPf^ATTA  TL  V PL  1.'^  (1  . -P  4 SF  APF4/pCF.  A P F A I / '■jrf  <;  F v^LIP'F 

ALPHAS  AN-GLF  FIF  ATTAC'<  - OEP- 
RRTPMsPPATTA IL  I IFT  FACTOR 
AL  PHAsO  ,r 

DATA  Vni.RA/  o.onp. 


VPLPAT 

PpTCV 

VPLPAT 

rhtpm 

VPLPAT 

ppTr> 

1 

0 , non , 

0.000  , 

O.?oo, 

-0.P15, 

o.aoo. 

-O.OFI  , 

2 

0 ,feno , 

-0.973, 

0 .Poc  , 

-1 .000, 

1.000, 

-1  . 000  , 

ALPHAsB. 

0 

3 

5.000  , 

vplrat 

RPTQM 

VPLPAT 

PPTfP 

VPLPAT 

pu  TP^ 

AJ 

0.000  , 

0.000  , 

0.200, 

- 0 . B 0 0 , 

0 , aoo , 

-0 . sa5. 

5 

O.feOO, 

-0.95P, 

0 .poo, 

-0 .opo  , 

1.000, 

-0.OP5, 

Al PHAsl  0 

.0 

6 

1 0.000  , 

VPLPAT 

PpTOKi 

VPLPAT 

PHTTU 

V P !,  P A T 

PPTO'K' 

7 

0.000  , 

0 . 0 0 0 , 

0 . ? 0 0 , 

— 0 . (1  h 5 , 

0 . u 0 (.  , 

-0..fiO?, 

8 

0.800, 

-0 . 90F , 

^.Mor , 

-o.ppp. 

1 . 0 0 0 , 

-0.930, 

ALPHAslP 

.0 

9 

15.000, 

VPLPAT 

PpTON 

VPLPAT 

PP  Tf '*>' 

V n.  P A T 

PPTP  H 

A 

0.000  , 

0 . 0 0 0 , 

O.2O0 , 

— 0 . /J  0 0 , 

0 , aO  0 , 

-0.69B, 

8 

0,600  , 

-0.«10, 

0 .Pon  , 

-0.«35, 

1 .ooo. 

-o.oao. 

Al PHAsPO 

.0 

C 

?0.0O0, 

VPLPAT 

RPTP^" 

VPL  RAT 

RP  Tt'*'' 

V r 1 P A T 

FPTP'" 

0 

0 . onn , 

0 . ooo , 

0,200, 

-0  , ?P0  , 

0 . a 0 r , 

-O.BO?, 

E 

0.6or,, 

-0.‘'30, 

0 . « 0 C , 

-0.671  , 

1 . 000  , 

-O.BPO, 

A 1 P H A s 3 0 

.0 

F 

30.000  , 

VPLPAT 

PPTPM 

VPLPAT 

PHTPF' 

VPI  PAT 

RpTr\i 

G 

0 . 000  , 

0 . n 0 0 , 

O.2O0, 

o.ooo. 

o.aoo. 

0,000, 

H 

0 .HOO  , 

n . ooo  , 

0 .600  , 

0.000, 

1 . 000  , 

0 . oco/ 
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siiPitf-sru  ir  ITU  r"PV(-  slppp  rP'^p-^vl  of  p? 


RFF.  Fir.,d.?.1  PFF.l 

BFMrRF  T 6 /KirSF  L/O 
RPFKrl  .n/RF^' 

F*FT,'srYl  IF'DFP  (L/r')/>.nsE  (L/OT 
ClArLTFT  c'l^VF  «:l'^FF  -d/Pan) 
F4FV  = C 

DAT#  rK'A2?l /n  . or  r , 


pFf 

Cl- A 

BF  ^ 

CL# 

PFN 

CLA 

2 

r .3r.r , 

1 .“30, 

0.900, 

1.900, 

0.600, 

1 .880 

3 

o.7rr , 

1 .P70, 

C .800  , 

1.860, 

O.ooo, 

1 .880 

a 

1 .oro. 

1 .880, 

F AF*'  rC 

.8 

s 

r.Rrn, 

«Fm 

Cl  A 

BFn, 

CL  A 

8F  N 

CLA 

h 

r.3rn, 

2.S80, 

0 . 0 0 , 

2.800, 

0.6O0  , 

2.620 

7 

n,7rn. 

2.810, 

0 ,8op , 

2.608, 

0.900, 

2. 800 

8 

1 ,oro. 

2.880, 

FAF\'=1 

.0 

9 

1 .nrr  . 

HFK. 

CLA 

RFN 

CLA 

BFF. 

CLA 

A 

r .300, 

2.890  , 

o.uoo. 

2.790, 

0.600, 

2.010 

8 

0.700, 

2.990, 

0.800, 

2.980, 

0.900, 

2.990 

C 

1.000, 

3.O00, 

FAF^ -p 

.0 

D 

2.000  , 

PF 

Cl  A 

RFr 

ri  A 

BFK 

ri.  A 

F 

0.300, 

2.730, 

C.uoo, 

2.860, 

0.600, 

3. 060 

F 

0.700, 

3.130, 

0.800, 

3. 200, 

0.900, 

3.290 

G 

l.Coo, 

3.290/ 

F A F N = 3 

.0 

DATA 

^^A232/3.0OO, 

PFK' 

CLA 

PF  f. 

OLA 

p F 

OLA 

1 

0.3O0, 

2.780, 

o.ooo , 

2.910, 

0.6O0, 

5.190 

2 

0.700  , 

3.220, 

0.800, 

3.290, 

0.900, 

3.360 

3 

1 . n 0 0 , 

3.^130, 

F AFM  = U 

.0 

a 

U.OOO  , 

pF  K 

Cl.  A 

PFK 

CLA 

RF  N 

rl.  A 

s 

0.300, 

2.770, 

0.900, 

2.920, 

0.800  , 

3.180 

h 

0.700  , 

3.290, 

0.800, 

3.3?o, 

0 . 9 r 0 , 

3.180 

7 l.nno,  ■^.ii7s, 

F J F N.  r R . r 

p B.ro(,, 
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TABLE  XIX  (Cont'd.  ) 


r 


BF^ 

ri.A 

RFN 

OLA 

MF  fv 

OLA 

0,300, 

?.7PO, 

0 . a 0 0 , 

?.«>ao. 

0 . F 0 0 , 

3.185 

0.700, 

3.?aF, 

O.POO, 

3.3?5, 

O.POO, 

3.ao8 

1.000. 

^.a75/ 

DATA  CNA2U/  0.000, 

0.300, 

1.850, 

o.POO, 

1 .P55, 

0,800 

1 .flF>0, 

o.fiOO, 

1 .880, 

1 . OOO, 

1.880, 

0,500, 

0,300, 

?.230, 

o.soo. 

? . a ? 0 , 

O.AOO 

0 ,000  , 

?.5ao. 

1 . 0 0 0 , 

2.580, 

1 . 000  , 

0.30O, 

?.580, 

0.500, 

2.7P0, 

n , POO 

?.«so. 

0 , BO  0 , 

2.‘)ao, 

1 .ooo. 

2.5P0, 

?,ooo. 

0.300, 

2. ‘>70, 

0.500, 

3.220, 

o.POO 

3.2«n, 

O.BOO, 

3.3?0, 

1 . ooo , 

3.200, 

3,000, 

O.300, 

3.230, 

0,500, 

3.a4JO, 

o.POO 

3.UQ0, 

0 .POO  , 

5.aoc, 

1 , 000  , 

3.a30, 

ii,  000  , 

O.300, 

3,3ao, 

0.500, 

3.530, 

O.POO 

3.':P0, 

o.POO, 

3.5P0 , 

1 . 000  , 

3.aF0, 

5,000  , 

0.300, 

.3.a35, 

O.50O, 

3.820, 

n , po  0 

3.^33, 

O.POO, 

3 . 8 0 c,  , 

1 . ooo , 

3.875/ 

TABLE  XX 

(sppf.  BcpK.  If  IJf-T  SL  f I'F  FHt-  fir,  I VF -C  VI  T Nf'E  RS 


FFP.  FTG.if.P.I  FFF.l 

RFMePFT  A/FiOSF  L/n 
P.PF^'sl  ,r/BF^' 

FAF^s  CVLINOpD  n/ni/MGSF  CL/^) 
rL4  = lTFT  riiPVF  Pi_npF-(’t/F/ir'i 

F A F f'  = 0 . r 


DATA 

OF.Apll/O  .000, 

RF  5 

CL  A 

PF» 

ri  A 

BF5 

CL  A 

1 

2 

0.200  , 

2.830, 

0.300  , 

2.810, 

0.800, 

2.7F0, 

0.500, 

2.380  , 

0,800, 

2. 380, 

O.POO, 

2,200, 

3 

1.000, 

2.220, 

F AFKrO 

.5 

u 

0.500  , 

OF  ^l 

OL  A 

HF‘' 

Cl  A 

RF  F 

CL  A 

5 

0.200  , 

2.780, 

0.300  , 

2.R1 0 , 

0.800  , 

2.P2^, 

8 

0.500  , 

2.835, 

0.800, 

2 . P 3 0 , 

O.POO, 

2. POO, 

7 

1 . 000  , 

2.780, 

F A F F.  = 1 

n 

• 

8 

1 ,O00, 

f F >• 

Cl  / 

RF  f 

ri  1 

• F > 

r 1 A 

9 

0.200, 

P.000, 

r . F 0 f, , 

2.020, 

0 . R , 

2.OF0  , 

A 

0.500, 

2.000, 

0 . p 0 0 , 

7.020, 

O.POO, 

7 . 0 fi  0 , 

B 

1 . 0 0 0 , 

F . oa5/ 
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Vtrrj  ^ 

DOi  /. 
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/ I . - -*2—  TABLE  XX  (Continued) 


DATA  CNA21 


RF,'- 

TL  A 

RFK 

01  A 

RFr. 

CLA 

(3  .?PO  , 

?.«PO, 

0.3O0, 

P.050, 

0 . u 0 ri , 

3.o?o 

o.soo. 

3,0<}0, 

0 . POO . 

3.iao, 

O.POO, 

3,P?o 

1.000, 

3.?Q0. 

FAF^'s3, 

0 

3.000  , 

RF*- 

n A 

PF  F 

CLA 

FFL 

Cl  A 

o.?oo. 

P.FA'O  , 

0.300, 

?.9R0, 

0 . R 0 0 , 

3.050 

r. , R 0 0 , 

3. 1 ?0, 

O.feOO, 

3.100, 

0 .ROO  , 

3.3P0 

1 . 000  , 

■«.390, 

F AF  F'Sii  , 

0 

U . 000  , 

RF^, 

CL  A 

PFN 

CLA 

pfn 

Cl  A 

o.?oo. 

?.PR0, 

0.3O0, 

P.930, 

0 , R 0 0 , 

3.030 

o.soo. 

3. 1 uo. 

O.POO, 

3.PP0, 

0 .POC , 

3.350 

1 . 0 0 0 , 

3. RIO/ 

FAFNrO.r 


PATA  C^A??/  P.OPn, 


PPF  Fi 

C!  A 

PPFR 

CL  A 

PPFF! 

CLA 

O.POO, 

1 ,P50, 

0.300, 

1 .530, 

0 . UOO , 

1 .730, 

0 . P 0 G , 
FAFRsO. 

O.500, 

1 .960, 

5 

O.POO, 

P.110, 

1.000, 

p.ppo. 

PPFK 

Cl  A 

PPFF 

Cl  a 

PPFF 

CLA 

0 . ?0  0 , 

1 .560, 

0 . 3 0 C , 

1 .900, 

0 . R p 0 , 

P.  1 po. 

0 . P 0 0 , 

F AF^.sl  . 
1.000, 

P.R70, 

O.POO, 

P.650, 

1.000, 

P.7bO, 

RPF^ 

CL  A 

PPFF 

CLA 

PPFN 

CLA 

0.200, 

1.630, 

C.300, 

P.ORO, 

0 .ROO  , 

?.3eo. 

0 .POO  , 
FAFm=2. 
2.0OO, 

P.P20, 

0 

0.800, 

3.020, 

1 .000, 

3 , ObO , 

PPF^ 

Cl  A 

BPFF' 

CLA 

PPFF 

CLA 

o.Poo, 

1 .7?o, 

0.300, 

P.200, 

O.ROO, 

P.PPO, 

0.600  , 
FAFMS3. 
3 . 0 0 0 , 

3.070, 

0 

O.POO, 

3.P90, 

1 .000, 

•^.?90, 

PPFK' 

CL  A 

PPFFi 

CL  A 

PPFF 

Cl  A 

o.poo. 

1 .790, 

0.300, 

3.2P0, 

O.ROO, 

P.F70, 

0 .poo , 

F AF  ^‘  = iJ  , 
U , 0 0 0 , 

3.  1 «o, 

0 

O.POO, 

3.300, 

1 . 0 r 0 , 

3.390  . 

PPFF 

Cl  A 

RPF  A 

01  A 

R P F 

rt.A 

O.2O0  , 

1 .«0O, 

0.300, 

P . 3 U 0 , 

O.ROO, 

P.FOO, 

O.POO, 

3.190, 

O'  . P 0 0 , 

7.R30, 

1.000, 

3 . R 1 r / 
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table  XXI 


m Ilf  TPn  fivF  IJFT  C'lt'VF-  SI  rpf- 

FPisAni'41  KinsF  FT^F^'FS.s  PiTin 
FC4Psf'nsF  E<|II^T‘■ESS  R4TT0 
Cl.  ARsFLI'^-f EO  CHMF  LTFT  curve  ‘IL'‘''f'E 
, FRlsJ .0 

DATA  CNA7P1/1 ,0fl0. 


FC  AP 

r q 

FCAP 

CfA« 

FCAP 

C v A R 

1 

O.ocn, 

1 .600  , 

o.ioc. 

1.670, 

0.200, 

1 .685, 

2 

P.3CC, 

1 .68.0, 

0.800, 

1,680, 

O.5O0, 

1 .590, 

3 

O.feCO, 

1 .'530  , 

0,700, 

1.850, 

0.800, 

1.380, 

ER1=1 .5 

a 

1 .5CP, 

FC  AP 

r^AP 

FCAP 

r I'AP 

FCAP 

CUAP 

5 

0.000, 

1 .800  , 

0.100, 

1 .875, 

0,200, 

1 , 865, 

b 

0,300, 

1 .790, 

0.800, 

1.785, 

0.500, 

1 .685, 

7 

0.600, 

1,608, 

0.7O0, 

1.850, 

0.800, 

1.^56, 

FRl s?.0 

8 

2.000, 

FC  AP 

CfiAR 

FCAP 

CNA« 

FCAP 

Of  AP 

9 

0.000, 

1 .885, 

0 . 1 0 0 , 

1 .887, 

O.2O0, 

1 .875, 

A 

0.300, 

1.8^3, 

0,800, 

1.795, 

0.500, 

1.718, 

B 

O.feOO, 

1 .6?0  , 

0.700, 

1 .500, 

0.800, 

1.365, 

FRl s3.0 

C 

3.O00, 

FC  AP 

CMP 

FCAP 

r A fi 

FC  AP 

rf:AR 

D 

0.000, 

1.9iJ5, 

o.loo. 

1.983, 

0.200, 

1 .922, 

E 

0.300, 

1 .882, 

0.8O0, 

1.818. 

0.500, 

1 .738, 

F 

0.600, 

1 .635, 

0.700, 

1.505, 

O.SOC, 

1 .358/ 

FPleii.O 

data 

CNA722/«.000 . 

FCAP 

CNAR 

FCAP 

r N.  4 R 

FO  AP 

r ^AH 

1 

0.000, 

1 .565, 

O.loo, 

1 .562, 

0.200, 

1 .980, 

2 

0.300, 

1.893, 

0.800, 

1.831, 

0.500, 

1 ,780  , 

3 

0.600, 

1.638, 

0 . 7 0 0 , 

1 .508, 

0,800, 

1 .358, 

FFI sR.O 

a 

B.ooo, 

fcap 

C^AP 

FCAP 

CMAP 

FO  AP 

ONAR 

5 

0.000, 

1 .980  , 

O.loo, 

1.972, 

0.200, 

1 .985, 

6 

0.300, 

1 .503, 

0 . 8 0 0 , 

1.835, 

0.500, 

1 .786, 

7 

0.600, 

1 .6ao , 

0.700, 

1 .509, 

0.8O0, 

1 .358, 

F01 S6.0 

8 

6.000, 

FCAP 

Cf.AM 

FCAP 

r fAR 

FCAP 

CNA8 

9 

0,000, 

1 .080  , 

0 . 1 0 0 , 

1 .975, 

0.200. 

1 .950, 

A 

0 . .3  I'  0 , 

1 , 5|I6  , 

0 . u 0 0 , 

1.838, 

0 . 5 0 0 , 

1 .788  , 

R 

0,600, 

1 .6U1  , 

0.700, 

1.510, 

0.800, 

1 .358, 

Fcir7.0 

C 

7 . 0 0 , 

F C AP 

C ‘ fi  P 

FCAP 

0 ' A R 

FO  AP 

r K A M 

D 

0 , 000  , 

1 .OOO  , 

^ . 1 0 .0  , 

1.978, 

1" . ? c , 

1 . 9501, 

F 

0.3O0. 

1.90'^, 

0 .g^O, 

J 

a 

• 

r .600, 

1.788, 

F 

o.Ri'o  , 

1 .6U?, 

■'.7^0., 

M 1 , 

r. . t.  n r, , 

' .168/ 
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TABLE  XXII 

HI  ijNTpn  nr,TVe  l.Tf-T  C'lPVF  SL  PPp 


ri  TA 


1 

? 

3 

41 


PA  T4 


8 

9 

A 

R 

C 

n 

F 

F 


FRlrATTlMl  ‘'CiFE  FI^F^'F^^8  PATjn 
FCAPs^pcF  PATTC 


r 1 A B = p 1. 1 1 ‘1 T F r 

nr,T  VF 

1 IFT  Cl'FVF 

SI  I'PF 

FFJSI .0 

rN<A731  f\  .OC'i, 

FT  A p 

r ► 1 B 

FC  A 1. 

OKA  H 

FOAP 

r M A 0 

0 . n n 0 , 

I 

.froo , 

0.100, 

1 . F 0 0 , 

0.200, 

1 .'^9. 

p.3no. 

1 

.FRO, 

r. . 41 0 0 , 

1 .SfrO, 

0.500, 

1 .530, 

r .8Pn , 

FPi =1  .F 

1 .SPO, 

1 

0.7O0, 

1 .ais. 

0.800. 

1.320, 

Fr  AP 

c^.  A« 

FCAP 

CNAP 

FOAP 

OMAR 

0 . OPO  . 

1 

.7fr?, 

0,100, 

1 .7fe1  , 

0.2O0, 

1 .756, 

P.3rn, 

1 

.7410, 

0,4100, 

1.709, 

O.5O0, 

1 ,650, 

0 . fr  0 rj , 

FR1 s?,n 
2.000, 

1 

.Fifr3, 

0.700, 

1 .4162, 

0.800, 

1 .335, 

FT  AP 

OFiAP 

FCAF 

C A R 

FOAP 

r 1-  A R 

0,000, 

1 

,«80, 

0.100, 

1.850, 

0.200, 

1 .8ai  , 

0.300, 

1 

.«1  ?, 

0 . 41  0 0 , 

1 ,7fr?. 

0.500, 

1 .693, 

O.frOO, 

FPIsT.O 

3.000, 

1 

,frOF, 

0.700, 

1 .a39. 

O.800, 

1.3a?, 

FT  AP 

r AiAP 

FCAP 

Cf^AR 

FCAP 

ri  AR 

0.000  , 

1 

.930  , 

0.100, 

1 .9?0, 

o.poo. 

1 .90F, 

0.3O0, 

1 

.nfr3, 

C,4J00, 

1 .ROP, 

0.500, 

1 .730, 

O.frOO, 

1 

.fr30. 

O.7O0, 

1 .50R, 

0.800, 

1 ,353/ 

FRIsU.O 

CNA732/U.OOO, 

F r AP 

^^'AP 

FCAF 

r MAR 

FOAP 

OMAB 

0.000, 

1 

.960, 

0,100, 

1 .950, 

0.200, 

1.930, 

0.300, 

1 

.ARP, 

O.aoo, 

1 .82?, 

.500  , 

1 .738, 

O.frOO, 

FPJ  sF. 0 
S . 0 0 0 , 

1 

.fr3P  , 

0.70C, 

1.51?, 

0,800, 

1 .3frO, 

FOAP 

r^■AP 

Ff  AP 

r >,AH 

FOAP 

r K AR 

0.000, 

1 

.97?, 

O.loo, 

1 .965, 

0.2O0, 

1 .9a?, 

0.300, 

1 

,900  , 

0,4100, 

1 .830, 

0.5OO, 

1 .7U5, 

O.frOO, 

FPi=H.n 

fr  . 0 0 0 , 

1 

. fr  0 0 , 

0.700, 

1 .'il  ?, 

0.800, 

1 .3fr?, 

FT  AP 

C^.AP 

FCAF 

TMAR 

FCAP 

'■MAR 

0.000, 

1 

.9Pfr, 

<".190, 

1 .975, 

O.2O0. 

1 .980, 

0.300, 

1 

.90?, 

O.aoo, 

1 .835, 

0.500, 

1 .7F0, 

O.frOO, 

FP1S7.0 

7.000, 

1 

.fr4l3. 

0.700, 

1 .51F, 

0,800, 

1.360, 

Fr  AP 

r » . A p 

FCAP 

OMAU 

FOAP 

AR 

0.000, 

1 

.9P0, 

0 . 1 OQ  , 

1.078, 

0.200, 

1 .-IFO, 

o.3rr , 

1 

.000  , 

0 . a 0 0 , 

1 .®3F. 

0 . F 0 0 , 

1 .7Fr , 

O.frOO, 

1 

'•■.70  0, 

O.frOO, 

1 . 
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Center  of  pressure  of  the  nose  alone  is  assumed  to  be  at  the  centroid  of 

the  nose  volume.  A correction  for  afterbody  ( ) is  based  on 

CPc 

data  of  References  9 and  10. 


) is  based  on 


Supersonic 

Nose  center  of  pressure  at  supersonic  speeds  is  based  on  data 
of  Figures  4.  2.  2.  l-18a  and  b of  Reference  1.  Table  XXIII  presents  center  of 
pressure  data  for  cone-cylinder  bodies  as  a function  of  Mach  number,  nose 
fineness  ratio,  and  the  ratio  of  afterbody  to  forebody  fineness  ratio. 

Table  XXIV  presents  similar  data  for  ogive-cylinder  or  bodies. 

Center  of  pressure  of  blunted  cones  and  ogives  is  based  on  data 
of  Figures  4.  2.  2.  1-21  and  4.  2.  2.  1-22  of  Reference  1.  Tables  XXV  and  XXVI 
present  center  of  pressure  for  blxmted  cones  and  ogives  as  a function  of 
nose  bluntness  ratio  and  actual  non-fineness  ratio.  At  Mach  numbers  grea- 
ter than  4,  the  center  of  pressure  is  based  on  a curve  fit  of  data  obtained 
from  modified  Newtonian  theory. 


List  of  Symbols 


CDgrp 


^Fi 

CLo^ 


Description 

Base  drag  coefficient 

Boattail  drag  coefficient 

Friction  drag  coefficient 

Induced  drag  coefficient 

Lug  drag  coefficient 

Zero -lift  drag  coefficient 

Skin  friction  coefficient  (corrected 
for  compressibility) 

Compressibility  correction  factor 

Incompressible  skin  friction  coefficient 
Lift  curve  slope 
Center  of  pressure 


A-39 


1 

\ 

V 

Des  cription  ; 

Nose  center  of  pressure  correction 
factor  to  account  for  afterbody  length 

Munk's  efficiency  factor  : 

Afterbody  (body  minus  nose)  fineness  ratio  : 

Launch  lug  equivalent  flat  plate  area  , 

Nose  fineness  ratio 

Body  fineness  ratio  { 

Boattail  lift  correlation  factor  based 
on  boattail  volume  and  geometry 

Boattail  length 

Cylindrical  center  section  length 
Missile  length 
Nose  length 
Mach  number 
Reference  area 
Body  wetted  area 
Volume  of  nose 
Angle  of  attack 
Beta  = M^  - 1 

5.  LIFTING  SURFACES 

The  lifting  surfaces  planform  may  be  either  the  basic  delta  or 
a clipped  delta.  The  fixed  lifting  surfaces  may  be  either  a planar  arrange- 
ment or  a cruciform  X arrangement.  The  movable  lifting  surfaces 
may  be  planar  with  a vertical  panel,  cruciform  +,  cruciform  X,  or  triform. 

Section  profile  may  be  double  wedge,  biconvex,  or  modified  double  wedge. 

5 . 1 Drag 

Drag  of  the  lifting  surfaces  consists  of  pressure  drag  due  to 
thickness  cind  skin  friction  drag 

Cq  = Cdp  -h 


J 


Symbol 

cPab 

CPq 

EK 

FA 

FL 

FN 

FRB 

^BT 


^BT 

LC 

LM 

LN 

M 

Ref 

Swet  B 
V 

a 
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TABLE  XXIII 

CPMP-r  YL  UOf- » fPK'TFB  PF  POFSS'lkF  ( 1C  ) 


BFMeBFTi/NOSF  FINFN'ESS  »ATin 
XCPDLsC.P.  IM  PFWCFNT  flnOY  LFMGTH 
FAF^*AFTEPB^^Y  (L/nV»''3SF  (L/H) 

I FAFMsI.fi 

DATA  XCPPP^/l .OCO, 

PFM  YCPDI  PF^o  XCPPI.  FF^  XCPDL 

1 0.?''0,  r\,j.07,  O.BOfl,  n.aFu,  o.uhT, 

2 1.000,  o,a7?, 

F AFK)s?,0 

3 2. OOP, 

PFN  XCPCL  pp^'  XCPOl  HFN  YCPOL 

0.?00,  0.?5«,  O.hon,  0,33?,  ('.BOC,  0,3?0, 

1.000,  0,3feB, 

PAFXia3.f' 

6 3.000, 

PFN  YCPPl  PP''  YCFPL  PFN  VCPOL 

7 0.200,  O.fcOO,  0.?6S,  C.BOO,  0,?BP, 

8 1 . 000  , i^.POT, 

FAFMsS.O 

9 5.000, 

PFN;  YCPOL  PP^  YCPHL  PP^J  XCPOL 

A 0.200  , O.lpix,  O.pno,  O.lPO,  0.800  , 0.1  97, 

8 l.OOO,  0.PO9/ 

F 4FM=1 , 0 

DATA  yr.pppp/i  , 000  , 

HPPFK  VCPP'  eoPFXi  YCPOC  PPFP^'  YPPPL 

1 O.poo,  0,aS3,  o.o'''',  0.478,  0,800,  O./IRO, 

2 1,000.  O.U7F, 

FAFM=?,o 

3 2.000, 

PPPFN  XCPOI  PPPFX  YOPPL  PPPFM  XOPPL 

a 0.200,  o,3B3,  0.800,  0.^88,  0.800,  0.37S, 

5 l.OOO,  0.38'', 

F AFX'x3 . '' 

8 3 . 0 0 0 , 

HFPFM  YCPPl  HPPFX  XOPDL  PPFFN  XrPOL 

7 O.2O0,  0.330,  0.8OO,  0.3)F,  0.800,  0.308, 

8 1,000,  0,?97, 


F4FK,ap.O 

9 5. cop, 

BPPF4  xcpni  opPFK  xrpC'i  ppfF''  xrp8[. 


( ' • 

u > 


r 
jr  t 


TABLE  XXIV 

nr.ivF  CYLTN.nfR  rPKiTFR  of  prfssupe  (supfr^i'mci 


FaF‘ si ,0 
FAF»>‘sl  ,0 

DATA  XCPFiJ/  1,000, 


PFN  YCPDL 

BFN 

XCPHL 

P,F  N 

1 

0.200,  0,?PiJ, 

0.900, 

0.335, 

0.700, 

2 

1.000,  0.3P0, 

FAF^■s?,o 

3 

2.000, 

PFN  VrODI 

BFF’ 

XOPPL 

PFn 

u 

0.?00,  0,i7fl, 

0 .900  , 

0.232, 

0,700, 

5 

1.000,  0.?8«, 

F AFN;S3,0 

6 

3.00O, 

PFM  XrODL 

fiFFi 

XCPOL 

PFNi 

7 

O.200,  o.ipF, 

0.900, 

0.17P, 

0.700, 

e 

l.ono,  o.?33. 

II 

• 

p 

U . 0 0 0 , 

FF\  XCPPI 

BFF 

xcphl. 

RF  N 

A 

0.20'',  O.OPS, 

0.900, 

0.199, 

0.7ro, 

M 

1.0  or,,  0.107, 

F AF  ^.5F  , 0 

C 

S.ooo, 

PFNi  XCPDL 

BFFj 

xcphl 

RFN 

D 

0,200,  O.OPO, 

0.900, 

0.121, 

0 .700 , 

F 

1 . 000  , n.i(S7/ 

FAF^sl .0 

DATA 

xcp(»f/  l.ono. 

RPPFn  YCPDI 

PPPF  Fl 

XCPDL 

ePPFM 

1 

0,300,  0.333, 

0 .500  , 

0.373, 

0.700, 

2 

1,000,  0,3FO. 

F AFMsP.O 

3 

2,000, 

PPPFm  XOPHL 

PPPF  N 

XCPDL 

PPPF  N 

a 

0.300,  0,?77, 

0.500  , 

0.295, 

0.700, 

5 

1.000,  O.PPS, 

FAF*.s3.0 

6 

3.0O0, 

BPPF^  XCPDL 

HPPFFJ 

XCPDL 

PPPF  K 

7 

0.300,  0.235, 

0.500, 

0.29F, 

0.700, 

F 

1,000,  0.232, 

F AFMsij  , 0 

9 

u.ooo. 

RFPFM  XCPDI- 

PPPF  N 

XCPDL 

PpPFFi 

A 

0.300,  0.200, 

0.500  , 

0.P1 3, 

0.700, 

B 

1,000,  0.197, 

FAFF's*.  , 0 

C 

5.000, 

ooPF^.  VCOp,  1 

PPPFM 

XOPIM, 

P P P F' ' 1 

0 

0.300,  0.177, 

0.500  , 

0. 1 PF, 

0 . 7 0 0 , 

F 

l.roo,  0,1 F7/ 

vr  pnt 

n,ihk. 


xrpHL 

0,?70. 


xcppl 

O.PlS, 


vrPOL 
n, 1 70, 


XCPOL 

0.150, 


yopi^L 

0.?P1  , 


xrPDL 

0,?07, 


XCMOL 

o.paT, 


xn  "OL 
0 , ?0«  , 


xrp^l 

0 , 1 FI  , 
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TABLEXXV 

RL'J^Ter.  c^'^F  rf-’TFt;  np  ppfvqi'pp  ( supf  B<5r^l  jc  i 


FCftPs  moSF  PL' SS  PATyn 
PPts  ArrUAL  (-ATTn 

XfP^sr.p,  c'P  ' P?r  I‘  p^oc?^T  ATTijiL  ' nsF  i ength 
1 = 1.0 

data  XCPNCl/l  .<^00, 


F r A p X r F 

FOAP 

XCP^ 

FCAF’ 

vCPf' 

1 

0,0''!^,  0 ,p:^o  . 

f'  . 1 0 f , 

O.P05, 

0.200  , 

A . 7P0 , 

2 

0.3''r,  f'.7PX. 

0 . a rn  , 

0.7?7, 

n .POO , 

A , FOB , 

3 

n.feoo,  ri,P7n, 

0.7^0, 

A ,^aA , 

A.BAO  , 

O.Fl  X, 

FPl=2.ft 

2.000, 

FTAP  XCP^- 

FCAP 

Xfpfj 

FCAP 

yf.PA' 

5 

0,000  , 

O.loO, 

0 . 8 a 2 , 

0 .200  , 

O.PIP, 

6 

0,300  , o.*;p9. 

o.aoo. 

O.PPO, 

o.pon. 

0.P2P, 

7 

o.pno,  o.aftS, 

n,7flo. 

0.a3B, 

0 ,P00  , 

0 . XBP  . 

FP1=3.0 

e 

3.0O0, 

FCAP  ycp^' 

FOAP 

XPPK' 

FOAP 

XCP^ 

9 

O.OOO,  0,PP7, 

r . 1 A 0 , 

0.P3P, 

o.poo. 

O.Bl  1 . 

A 

0.300,  o,5«n. 

0 .aofl , 

O.P'JP, 

0 . P 0 0 , 

0.P1 0, 

B 

O.feoo,  n.opp. 

0 . 7 0 0 , 

o.aio. 

A. BOO, 

o.X3», 

FR{sU,f) 

C 

U , 000  , 

FTAP  XCDM 

FCAP 

XCOM 

FCAP 

yf  Ph 

D 

0,000,  0,P6P, 

o.ion. 

0.‘.37, 

■ •>.2  0 0 . 

0 . B 1 A , 

F 

0.3'^'^,  n.P7p, 

0 .ann , 

0 .pax. 

0 . P 0 A , 

A . FOX , 

F 

O.POO,  O.UPP, 

0.7fl0, 

n.39P, 

o.«ro. 

O.Xi?1  / 

FP185.0 

DATA 

XCPMr2/5.000, 

FTAP  XCP^i 

FCAP 

X c 

FCAP 

yr  P'-i 

1 

0,000,  0, hhh , 

0,100, 

0.P3P, 

o.?-"0. 

O.Bfi7, 

2 

0.300,  O.P73, 

o.aon. 

O.P3Q, 

o.poo. 

A , a 9 B , 

3 

O.PPO,  o.aaP, 

r . 7 0 r , 

o,39a. 

0 . F 0 r , 

0 . XI  1 , 

FPtsP.O 

a 

6.0ft0, 

FOAP  xCPf 

FC/P 

yr  PKi 

FCAP 

yCPf'J 

5 

o.on'', 

0 . 1 0 A , 

o.P3‘=, 

0.200, 

0 .POA), 

o.3ro,  n.P7?, 

0 . a r 0 , 

O.S3P, 

o.poo. 

0 . a 9 p , 

7 

' 0 . 6 ''  0 , 0 , a /I  a , 

r . 7 A 0 , 

0 .3PP, 

O.FIAA, 

O.XAF, 

Fcis7,n 

P 

7 . 0 il  0 , 

Fra®  xcP’' 

F TA  P 

« r p * ' 

FT  A P 

Yf  P'l 

9 

r,  ,non,  'i.Ph*., 

0 , 1 0 r , 

0 . ? 0 A . 

0 . B I-  F . 

A 

0 .3P(‘  , f'  .P70  , 

r . a A r , 

0 .p  i'=;. 

A . R r A , 

' . a u X , 

R 

n.pro,  n.aax. 

. 7 A r, . 
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TABLE  XXVI 

R1  IlMTPn  r'r,l\/s:  t■F^TF°  of  F»p«;s,icr  f SllPPWenf  ir  I 


FFAPs  ‘-nsF  pi,ii‘Tnfss  PATin 

FP1=  /\rTiiAL  r.nsF  fjA'FAiESS  patK' 


DATA 

XCP^sr.P.  f’F  ^'HRF 
FU1=1 .D 
VCP7F)/I.nnn, 

IM  pftvf^p.T 

actual 

►nSF  1 F.fvGT.^ 

FTAF  yCPl" 

FCAP 

XCP'' 

FCAP 

yC  Pu 

o.oon,  n.uqis. 

P . 1 P P . 

O.Ufl7, 

p.pppj. 

P.U7F, 

ri^UhTi, 

R .UPP  , 

P.U51  , 

0.5PP, 

P .U3®  . 

O.Dfin,  ('.a??, 

FRlrP.R 

? . n p 0 . 

P .7PP  , 

P.UPF, 

0 . e p 0 , 

P.3P1  , 

F r A c y c P ^ 

FCAP 

xr  p*'! 

FCAP 

vCPn 

p.ppp,  o.«J7n, 

P.lOP, 

P.UR7, 

P.PPO, 

P , uu 1 , 

o.3np.  P.ii?E>, 

P.ilPfi, 

P.UP9, 

P.5PP, 

P.3RP, 

P .bpr  , P , 

F R 1 r 3 . P 

3 . P P P , 

P . 7 P 0 , 

P.31S, 

p.flpp. 

P,PR7, 

FCAP  yfP'. 

FCAP 

ycpv 

FCAP 

yCPi'* 

P.PPP,  P.RDO, 

P . 1 P P , 

P.U53. 

0 . 2 P P , 

P.U3‘', 

p.3rp,  P.RJP, 

P .UPP  , 

P . UPP  , 

P.SPP, 

P.37Q, 

P,6PP.  P.TST, 

F p 1 = n . p 
u.opn. 

P.7P0, 

P.3)7, 

P.PPP. 

P,?7T, 

Fr«p  ycpfi 

FCAP 

ypRu 

FCAP 

yr  p\ 

P.PPP.  P , 

O.IPP, 

P.US1 , 

P.PPP. 

P.U33, 

P.3PP,  P.UIP. 

P.UPO, 

P.3°u, 

P .FPP  , 

P,371  , 

DATA 

P.FPP,  P.3<JR. 

FRlsS.C 
yCP7S?/S.P0P. 

0.7PP, 

P.3PS, 

p . fiP  P , 

P . ?t.P/ 

ftap  ypPA' 

FCAP 

ypPM 

FCAP 

ypPA. 

P.PPP.  P.RF,7, 

P . 1 PP  , 

P.U5P, 

P.PPP, 

P.U30, 

P.3PP,  P.P1  p, 

P . 0 P C , 

p.3n9. 

P.‘=PP, 

P.3FU. 

P.6PP,  p,T3F, 

FPlrD.P 

A^.OPP, 

P .7PP  , 

P.PPP. 

P.fiPP, 

P.PFO, 

F r A p y c P 'J 

Ff  AP 

ypp^ 

FCAR 

y r P : . 

P.PPP,  p .4iFs7, 

P.IPP , 

P . UUR , 

P.PPP, 

P.U?", 

P.3PP,  p.aps. 

p.app. 

P.3UC, 

P.^PP, 

P.3^P, 

P.F>PP,  P.33P, 

F p 1 r 7 . P 

7 . C P P , 

P.7PP , 

P.PPI , 

P . F P P , 

P.?«7, 

Fcap  yppf 

FCAP 

yrP^i 

FC  A F 

VC  Pf" 

P.PPP.  P.UF.7, 

P.lPP, 

P . uu  7 , 

P.?PP. 

P.U?7, 

P.3PP,  P.RPR, 

P . UP  P , 

P.3R?, 

P.srp, 

P.3SA, 

P.PPP,  P.3?R, 

P.7PP, 

A-44 

p.P^P, 

p . ® P P , 

p . P3U/ 

5.1.1 


Zero-Lift  Drag  Coefficient 


Subsonic 

The  lifting  surface  drag  at  subsonic  speeds  is  based  on  skin 
friction  drag  corrected  for  thickness  ratio  effects.  The  skin  friction  coef- 
ficient is  based  on  Reynold's  number  computed  from  the  mean  geometric 
chord  of  the  surface.  Lifting  surface  zero  lift  drag  coefficient  at  subsonic 
speeds  is  given  by  the  following; 

Surface  wetted  area  is  assumed  to  be  twice  the  planform  area  S . 

P 

N accounts  for  the  number  of  panels. 

P 

Transonic,  M = 1 

Zero-lift  drag  in  the  transonic  region  consists  of  friction  drag 
plus  a pressure  drag  term  derived  from  a curve  fit  of  data  of  Reference  2. 


Supersonic 

Zero-lift  drag  coefficient  in  the  supersonic  region  consists  of 
friction  drag  and  pressure  drag  due  to  thickness.  For  subsonic  leading 
edges,  zero-lift  drag  coefficient  is  given  by  the  following: 

For  supersonic  leading  edges 

The  airfoil  section  constant,  K,  is  based  on  data  of  page  4.  1.5.  1-15  of 


Reference  1. 


5.1.2 


5.1.3 


Induced  Drag  Coefficient 


Induced  drag  is  based  on  the  relation 

Cp^  = ^ 

Total  Drag  Coefficient 


Total  drag  coefficient  consists  of  zero-lift  drag  coefficient  plus 
the  induced  drag  coefficient 
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5.  2.  Lift  Curve  Slopes 

Subsonic 

Isolated  lifting  surface  lift  curve  slope  is  obtained  from  data 
of  Reference  1 1 . 

The  arrangement  factor  Kp  accounts  for  the  reduced  effective  angle  of 
attack  and  rotates  the  lift  vector  into  the  pitch  plume  for  tri-form  and 
cruciform  "X"  configurations.  Table  XXVII  presents  isolated  surface  lift  curve 
slope  divided  by  aspect  ratio  ^Yoc//^  ^ function  of  surface  taper  ratio  (A.)i 

the  product  of  aspect  ratio  eind  the  tangent  of  the  midchord  sweep  angle 
A and  the  product  of  aspect  ratio  and|S  \lM^  - 1.  (BAR) 

Transonic 

Lift  curve  slope  of  the  isolated  panel  is  defined  by  the  relation 

A.  , (3)^^ 

Supersonic 

Lift  curve  slope  of  lifting  surfaces  with  exposed  aspect  ratios 
greater  than  1. 0 is  based  on  the  relation. 

-/("A  ^ 

Low  aspect  ratio  wings  (A^l.  0)  produce  more  lift  than  indicated 

by  Reference  H according  to  data  of  Reference  10  and  12.  A correction 

factor  (1  + K K, ,)  was  developed  by  Martin  to  account  for  aspect  ratio  and 
AM 

Mach  number.  This  factor  is  applied  to  the  obtained  from  Reference  H. 

' |S)]  ( • ^ 

5 . 3 Interference  Effects 

The  mutual  interference  effects  of  the  wing/body,  tail/body,  and 
tail  in  the  presence  of  the  wing  are  based  on  the  methods  described  in  Ref- 
erences 1 and  13. 
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TABLE  XXVII 


SURSHmTC  SUBPaCF  I IFT  f|.cvF  Sl  4 r F /s.spftt  patTp 


TPIsSIIPFacF  tapfr  patT' 


DATA  A ir^P  1 /n , non , no.?sn,  o.'^nr. 

1 ,009, 

ATKS?I  = 4?PFrT  MATTO  TIa-FS  TA^^.F^T 

nF  Pircw 

1 

o.noo, 

1.000, 

2.O0O,  3.00O, 

. 0 n 0 , 

5. 900  , 

RA PIsOETA 

TI'-FS  A8PFCT  PATIO 

2 

0 , nno , 

0.500. 

1.O0O,  1.5O0, 

2.909, 

3 

3.000, 

a.  000  , 

F.OOO,  6.0O0, 

7.009, 

ATCL  IsSUPPirE  LIPT  C'JPvF  51. 

nPF/ASPFCT  patT 

TPJsO . 0 

4T‘:??I  = 0.0 

a 

1 .570, 

1.500, 

I.a?5,  1.3P0, 

1.219, 

5 

1.030, 

0.900, 

0.780,  n.700. 

0.669, 

TPTxO.O 

ATA.SPIsl  .0 

b 

1 .570, 

1 .<i90. 

1.^00,  1.285, 

1.170, 

7 

1.010, 

0.880, 

0.770,  C .690, 

0.629, 

TRTsO.O 

AT^8?I=2.0 

6 

1.570, 

1 .«oo. 

1.310,  1 . 18C, 

1.080, 

9 

0.950, 

O.8U0. 

0.750,  0.670, 

9.610, 

TRTsO  , 0 

4TMS2I=3.0 

A 

1 .260, 

1.190, 

1.120,  1.055, 

0.990, 

B 

0 ,B90  , 

0.800  , 

0.720,  0.65r, 

9.800, 

tpTsO,  0 

ATAIS2ISP.0 

C 

1.0U5, 

1.015, 

0.975,  0.935, 

0.899, 

0 

O.BIO, 

0.7U0, 

0.670,  0.620, 

9.8«0, 

TPtcO.O 

4TMS2I=5.0 

E 

0.900, 

0.888, 

0.880,  0.883, 

0.009, 

F 

o.7ao. 

0.600  , 

0.6a0,  0.590, 

0.550, 

TRIsO.O 

AT‘ S2I=6.0 

r. 

0.790, 

0.783, 

0.770,  0.750, 

0.730  , 

h 

0.680, 

0.630, 

0.590,  0.560, 

0.520/ 

TRT  = 0,?5  aU\S2T 

= 0.0 

DATA  Air602/1 .570, 

1.520,  i.aao. 

1 .365, 

1.250. 

1 

1.070, 

o.93« , n.82C , 

0.730  , 

0.680, 

tpt=o.?5  AT^52I 

= 1.0 

2 

1 .570  , 

1 . 5 1 , 1 . u 2 r. , 

1.325, 

1.220, 

3 

1 .0£in  , 

0.P1  0,  r .791^, 

0.720  , 

0,650, 

TPT  = 0.28  ATK  SPI 

= 2.9 

u 

1 .570  , 

1 .PU9  , 1 . 330  , 

1 .229  , 

1 .1  ■»5, 

5 

O.PfiO, 

0.87'^,  f , 7 7 , 

9.899, 

fi . 8 3 9 , 

TP  T = r , 25  A T ► 62  T 

= 3.0 

b 

1.265, 

1.209,  1 . 1 ^ 9 , 

1 .088, 

1 .990. 

7 

0.910, 

9.81'',  r . 730  , 

9.8^9, 

0 . 1 <■  , 

TPT=9,28  iTA?2I 

= 8.' 

1 .080,  , 

1.919,  A.qyr,  . 

9 . 9 3 9 , 

9.890. 

9 

n . 8 X 0 , 

9 .■'69  . 9 ^(.or  , 

0.869, 

9.80  - , 

p.onr. 


o o 


1 . 


TABLE  XXVII  (Contipurri) 


TPTsp 

.?5  4T^P?I=S.O 

4 

n . 9 f)  P , 

O.P70, 

0,P4O,  O.SJO, 

0.7P0, 

B 

0.740, 

0,890, 

0.840.  0.800, 

C.980, 

T » T s 0 

.?5  41^S?T=8,0 

c 

n.79P, 

0.784, 

0,740,  o.7?r. 

0 . 7 0 0 , 

D 

0.870, 

0.840, 

0 ,80r  , 0.480  , 

0.930/ 

TPTsO 

.40  4T9S?1=0.0 

04T4 

410803/1 .87P, 

1 .5?o . 

1.430,  1.340, 

t .?40, 

1 

1 . OPO , 

o.oto. 

n.8?0,  0.730, 

0,880, 

T C T = A 

.40  AT^5?I=1.0 

? 

1 .=■70. 

1 .900, 

1 .40S,  1 .700, 

1 . ? C 0 , 

3 

1 . 0 4 P , 

0.9?0. 

0.800,  0.7?n, 

0,840. 

TOT  = n 

,40  4 T^:S?I=<J,0 

4 

1 .=■70, 

1 .430  , 

1 .31  Of  1 .?00. 

1 .1?0. 

S 

0.940, 

0.B70, 

0,770,  O.^qo, 

0.870, 

TPIsO 

.90  4T^'S?I=3.0 

t> 

1 .?P0, 

1 .?00 , 

1,130.  1 .080, 

1.010, 

7 

0.910. 

0,p?0, 

0.730,  0.880, 

0.800. 

T P T 5 0 

,40  4TM4?I=4,0 

ft 

1 .080, 

1,017, 

0.979,  0.940, 

0.900, 

9 

o.p-^o. 

0.7S0, 

0,8PC,  0.8?0, 

0.5P0, 

TPTsO 

.40  4TNS?T=5, 0 

4 

0.91 0, 

O.P74. 

r.H40,  0.P14, 

0.7P9, 

0.740, 

0.8P0  . 

0.830,  0.990, 

0.990, 

rPJzP 

.4P  4 7>'5?J  = 8.0 

r 

0.7n0, 

0,780. 

0.740,  0.719, 

0.890, 

n 

0,880, 

0.830, 

0.980,  0.990, 

r .9?o/ 

T 4 T = 1 

,op  4T''S?I  = 0,O 

04  7 4 

4 10804/  1 .970, 

1 .9?=. 

1.449,  1.340, 

1 .2?0, 

1 

1.040, 

0 . 90  0 , 

0.790,  0.700, 

O.840, 

T9  Ts  1 

,oc  4T^5?I=1.0 

? 

1 .S70, 

1.910, 

1.4?0,  1.309. 

1.190, 

3 

1 . 0 0 0 , 

0 , apo , 

0.770,  0.890, 

0 , 8 ? 0 , 

TPIs1 

.00  ATK.4?I=?,0 

u 

1.570, 

1.477, 

1.390,  1.??0, 

1.108, 

5 

0.940, 

O.PPO, 

0.730,  0.880, 

0.810, 

TPjsl 

.oc  A7''P?I  = 3.o 

h 

1 ,2P0, 

1.170, 

1.093,  1.o?n, 

0.980, 

7 

0 , P8  0 , 

0.770, 

O.890,  0.830, 

0.9P0, 

T p T = 1 

.00  4T8P?I=4,0 

n 

1 .040, 

0,99?. 

o.9ur.  0.aR8, 

0,849, 

9 

0. 7«0, 

0.71 0, 

0,840,  0.800, 

0..98C, 

TP|  = 1 

.00  iTVS?I=9,0 

4 

0.9OS, 

o.«83. 

C .P03,  0.7,,8, 

'■'.  740  , 

B 

0.700, 

".840, 

0,800,  0.480, 

' .940, 

TP  T = 1 

.00  iTK,4?T=8,0 

c 

0,790, 

0.740, 

0.709,  0.870, 

0.840  , 

0 

0.8?0, 

0,990, 

0.980  . O.c-so, 

0 . 900  / 
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TABIjE  XXVII  (Continued) 

SUPfRSOK'lC  Sl'PFiCP  I TFT  niPVf-  SlC  PF /S'lpp  if.F  ASPFFT  PtTTn 


TP IsSiipfacf  taper  rattm 
data  AIC61 1/n.onn,  o.PSC,  o.snn,  i.ooo, 

AT\??TaASPFrT  RATir  TT^FS  Tant-emt  pF  >^irr,-r,pr)  S'^'FFP 

1 0.000,  1.000,  ?,000,  3. ''0  0,  u.onn, 

2 5.000,  6.000, 

RAPIsPFTa  TI^'FS  aspfct  RaTTH 


5 

0.000, 

0 

.500, 

1 . 0 0 c , 1 . 5 0 0 , 

2 

. 000  , 

u 

3.000, 

a 

. 000  , 

5.000,  6,000, 

7 

. 000  , 

AICI 1= 

SURFAOF  l 

TFT  rURVF  5LnPE/ASPF,r.T 

H A T T r, 

TRTs 

0 

,00  ATtJSPIsO.O 

5 

1 .570, 

1 

.700, 

1.750,  1,730, 

1 

.700, 

6 

1.210, 

0 

.950, 

0.770,  0.650, 

0 

.550  , 

TRTs 

0 

.00  ATMS?I=1.C 

7 

1.570, 

1 

.7ao, 

1.7a5,  1.630, 

1 

.a90. 

8 

1.250, 

0 

,960, 

0.790,  0,660, 

0 

. 560  . 

TRTs 

0 

,00  ATM.5?Ts?.0 

9 

1.570, 

1 

.530, 

l.a60,  1.160, 

1 

.300  , 

A 

l.UO, 

1 

.010, 

c.e?o,  o.iRHo, 

0 

.5«0, 

TRTs 

0 

.00  ATMS21=3,C 

P 

1 .260, 

1 

.300, 

1.275,  1.230, 

1 

.170, 

c 

1.050, 

0 

.950, 

0,660,  0.710, 

0 

.600  , 

TRTs 

0 

.00  AT  F' 821  = 9.0 

D 

1 ,0410, 

1 

.066, 

1.060,  1.075, 

1 

. 060  , 

F 

0.9«0, 

0 

.900  , 

0,630,  0.760, 

0 

.630  , 

TRTs 

0 

,00  ATRiSPTsS.C 

F 

0.900, 

0 

.515, 

0.930,  0.93R, 

0 

.oao, 

G 

0.910, 

0 

.flao, 

0.760,  0.7 30, 

0 

.68  0 , 

TRTs 

0 

,00  ATKSPIsh.O, 

H 

0,600, 

0 

.800  , 

O.6O0,  0,6  on. 

0 

. 800  , 

X 

0.600, 

0 

.610, 

0.750,  0.710. 

0 

. 660  / 

TRTs 

0 

.?5  AT^-SPfsO.O 

data 

A 1061  2/1  .570  , 

1 

. aoo , 

1.690,  1.630, 

1 

.680, 

1 

1 .2P0 , 

0 

.9/iO  , 

0.77O , 0.650  , 

0 

.560  , 

TRTs 

0 

,25  AT^'S2Isl.0 

2 

1.570, 

1 

.350  , 

1.6ao,  1.710, 

1 

.560  , 

3 

1.220, 

0 

,9ao. 

0.77C,  0.650, 

0 

.560  , 

TRTs 

(■ 

,?6  AT>-:6?J=?.0 

a 

1.570, 

1 

.120, 

l.a60,  I.a70, 

1 

. a 1 0 , 

5 

1.210, 

'1 

. oqn  , 

0.810,  0 . 8 0 . 

-6 

.500  , 

TRTs 

0 

.R6  AT>  FPT S3. 0 

h 

1.266. 

1 

.160, 

i.nQ(, , 1.280, 

1 

.26  ' , 

7 

1.130, 

1 

. 020  , 

^.680,  O,7e0, 

r, 

.51  ' , 

TRTs 

r, 

. ? 6 A T ‘ 8 P ] = 8 , 0 

A 

1 5V, 

.061, 

1 . f'  6 r , 1 . 0 R . 

1 

, r ii  , 

9 

1 . a 0 , 

0 

.9R0. 

'■.87r,  ''.16  0, 

r 

. ^ • 
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TABLE  XXVII  (Continued) 


TPTsn 

,25  ATK'<;2J=5.8 

A 

0 . non  , 

.°)  8, 

8.928,  8,^30, 

8 , <3  a 8 , 

H 

.oto  , 

8 

.0)8, 

8.P58,  8,7>18, 

0.670. 

TP  Ts8 

,26  ATMC2T=h.8 

c 

8 

.7PR, 

0.7n(;,  8,78  8, 

8.7P0, 

n 

'^.800, 

n 

.^38, 

8.P20,  8,768, 

C . 7 8 0 / 

TPT  = 8 

,R8  ATKiS2I=8.0 

DATA 

4Tr^i3/i.R7o. 

1 

. «|08  , 

1 .960,  1 ,R30, 

1 .590, 

1 

1 ,2f'0. 

n 

.038  , 

0,760,  8,f,U8, 

8.560, 

TWT  = 8 

,60  AT6S2Ts).8 

2 

1 .57^, 

1 

. 7Pn  , 

1.9  TO,  1.7B0, 

1 .‘*70. 

1 

1 .Pnn . 

8 

. OU  8 , 

8.778  , 0 .668  , 

8.560., 

TP  T = 8 

,60  ATMS?I=2.0 

a 

t .S70, 

1 

. a P 8 , 

1.360,  1.a38, 

1.958, 

S 

1 .21  'I, 

8 

.0«8, 

8,P0C,  0.670, 

0.5P8, 

TPTsO 

.50  A TK)S2I=3.0 

6 

1 .278, 

1 

1.2?5, 

1 .255, 

7 

1 . 1 '•r  , 

1 

.818, 

O.pao,  8,708, 

0.680, 

TPT  = 8 

,60  AT6CPT=0.C 

8 

1 . 8 ^ n . 

1 

. 868  , 

1,060,  1.868, 

1.868, 

0 

1 

8 

.0P8  , 

r.P78,  n.Tso, 

0.6/10  , 

TC  T = 8 

,50  A7^'5?Tr5,8 

A 

n.Rin, 

8 

.018, 

8.918,  8.018, 

8.910, 

R 

0 . « 2 8 , 

8 

.018, 

f,  .P68  , 8.778, 

^.670, 

TRT  = 8 

.50  ATMS2Ts6,0 

c 

8. 78  8, 

8 

.708, 

0.P80,  O.ai 8, 

O.B?0, 

D 

8 

.1*38, 

0 . P 2 8 , 8,778, 

0.780/ 

T«Trl 

,88  AT^P2IsO,0 

DATA 

A7C8ia/i.57n, 

1 

.798  , 

1.998,  1.«08, 

1 .‘=28, 

1 

1 . ) 28, 

8 

. PP  8 , 

8.738,  0 ,628 , 

8.5ao, 

TPTr  1 

. 8 c A T ti  s 2 I = 1 , 8 

? 

1 .'^70, 

1 

. 007  , 

2.018,  1.720, 

1 .aP8 , 

3 

1.138, 

8 

. 008  , 

0.7ao,  0.638, 

0 . 5 0 C , 

TPTr1 

,80  ATK'P2I=2.0 

a 

1 .578  , 

1 

.537, 

1 .528  , l./JPO, 

1 . a a 8 , 

s 

1.188, 

8 

.038  , 

8.768  , o.#.58. 

^.560, 

ThT  = 1 

.80  47652153,0 

8 

1 .?88, 

1 

. 388  , 

1.300,  1,275, 

1,290, 

7 

1.108, 

8 

.068, 

O.POO,  8,6®8, 

0.5BO, 

T P 7 = 1 

, 80  4T'^S2l=a,8 

8 

1 .0'=8, 

1 

. 858  , 

1 . (1 5 0 , 1 . 8 a 8 , 

1.820, 

9 

8.088, 

8 

. 078  , 

0.P30,  8,710, 

8.610, 

T»  T=  1 

.80  ATS  52  1=5,0 

A 

8.018, 

8 

, 088  , 

8 ,PB8,  8,B73, 

0 .«65. 

P 

8 , 8 ^ r.  , 

8 

. R68  , 

O.Pac,  8,7ao, 

0.668, 

TPT  = 1 

,08  4 7^.521=6,0 

C 

8.708, 

8 

. TMP  , 

0.7B8,  0,778, 

''.  765  . 

0 

8.788, 

8 

.758  , 

0.760,  8,7B8, 

8.728/ 
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5.3.1 


Carryover  IXie  to  An^le  of  Attack 

-f ) - ^ J - i-?3  ^ ^ f Jl 


e-fj 


^acs) 


The  basic  carryover  factors,  K and  K 

S(B)  B 


are  based  on 

slender  body  theory  as  presented  in  References  1 and  13.  At  supersonic 

speeds  the  carryover  from  the  surface  to  the  body  has  been  modified  to 

account  for  the  aft  movement  of  the  area  influenced  by  the  surface. 

5.3.2  Carryover  Due  to  Incidence 

The  carryover  factor  case  of  surface  deflection  is 

based  on  slender  body  theory  as  derived  in  Reference  13.  The  following 

expression  gives-^3  in  terms  of  the  semi-span-ratios  ratio  {T  = s/r) 

(B) 


J 

.._L 

-i-  TT  crSiV 

= 

r’^CTi-i) 

(r\,r 

T^(r-0 

L / •/  Z ^ 

[^,n  ' r -! 

I 

The  carryover 

factor^D 

is  based  on  i 

fr 


tion  that  the  surface  transmits  a certain  fraction  of  its  lift  to  the  body  whether 
the  lift  is  developed  by  angle  of  attack  or  incidence  angle.  An  approximate 
value  is  given  by  multiplying  the  value  of 


'(S) 

carryover  factor  due  to  angle  of  attack  Kg 


,/Kc 

(S)  S(B) 


c by  the  ratio  of  the 


5.  3.  3 


Asta) 

Wing  Vortices 

The  wing-tail  interference  results  from  downwash  in  the  region 
of  the  tail  caused  by  wing  vortices.  Figure  7 illustrates  the  vortex 
model  used  in  determination  of  the  wing -tail  interference.  The  lateral 
location  of  the  wing  vortices  is  based  on  slender  body  theory  and  is  given  by 
the  relation: 
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[ 


I 

i 


i 

f 


ri-rs-i-j/fC'-f/J-f 

Z (l-f 

For  low  angles  of  attack  and  zero  deflection  of  the  forward  surface,  the 
vertical  location  of  the  wing  trailing  vortex  is  assumed  to  lie  near  the  plane 
of  the  wing.  The  tail  interference  factor  "i  " is  evaluated  in  terms  of  the 
two  external  and  two  internal  or  "image"  vortices  as  defined  in  Reference 
14.  Appendix  B of  Reference  13  describes  the  use  of  strip  theory  to  obtain 
the  tail  interference  factor  from  integration  of  the  lift  contribution  of  the 
vortex  system  on  the  tail  panels.  The  lift  on  the  tail  section  due  to  wing 
vortices  is  computed  from  the  relation 

~ 

5. 4 Center  of  Pressure 
All  Mach  Numbers 

Lifting  surface  center  of  pressure  at  pf  = 0 is  obtained  from 
data  of  Reference  11. 

CP  = f (A,  X , tan  ^ > 

Table  XXVIII  presents  lifting  surface  center  of  pressure  data  in  percent  of 
mean  geometric  chord  as  a function  of  aspect  ratio,  taper  ratio,  tangent  of 
the  midchord  sweep. 

5 . 5 List  of  Symbols 


Symbol 

A 

b 


Description 
Aspect  ratio,  exposed 
Surface  span,  2 panels,  no  body 
Drag  coefficient 
Induced  drag  coefficient 
Zero-lift  drag  coefficient 
Skin  friction  coefficient 
Pressure  drag  coefficient 
Lift  curve  slope 
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Symbol 


Des  cription 


CP 

f 


Center  of  pressure 

Lateral  location  of  wing  trailing  vortex 


Tail  interference  factor 


Airfoil  section  constant 

Correction  factor  for  C^^  of  low  aspect 
ratio  surfaces  “*■ 


^B(S) 


Lifting  carryover  factor  due  to  angle  of 
attack,  surface  in  presence  of  body 


''b(S) 


Lifting  carryover  factor  due  to  surface 
incidence,  body  in  presence  of  surface 


Lifting  surface  arrangement  factor 


Correction  factor  for 
ratio  surfaces 


C of  low  aspect 


^S(B) 

''S(B) 


r 

s 


(t/ c) 


Lift  carryover  factor  due  to  angle  of  attack, 
surface  in  presence  of  body 

Lift  carryover  factor  due  to  surface  incidence, 
surface  in  presence  of  body 

body  radius 

Surface  semi-span 

Planform  area  of  1 exposed  panel 

Reference  area 

Thickness  ratio) 


Angle  of  attack 
- 1 ' 

Leading  edge  sweep  angle,  measured  from 
perpendicular  to  root  chord 


Mid-chord  sweep  angle,  measured  from 
perpendicular  to  root  chord 

Surface  taper  ratio 


r 


Semispan  - radius  ratio 
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TABLE  XXVIII 


ff-  f'.f.  LTt-TiMG  S'mcArF f si'Rsn^  ir ) 


TPIsTAFPW  CiTin 

DATA  yAr78i/0,P00,  O.PSO,  I.oon, 

ATM<;?T  = A<^Pf  TT  H&TTn  TI^FS  IAnGF^’T  [)P  Mirrv-ron  <?wPFP 

1 n.oor,  l.nr.o,  p.rnc,  3.nor. , u.non, 

hARIsASPEr.T  RATin  TI”PR  BFTA 

2 n.oor.,  o.sno,  l.non,  i.BOO,  ?.nor, 

3 3. nor,  a. 000,  B,0(i0,  is.ooc,  7. nor, 

XBATsBiiPF  Af  P rp  J *.  ppPCFNT  AN  ppOMPyPTf^  pufiPf) 
TRTsO.O  AT^,S?IrO.O 

B n.?F,o,  o.?B9.  O.PftB,  0.2FP,  0.27P, 

5 ').27li,  0,?73,  0.?71,  0,270,  0.?#!7, 

TPTsO.O  ATnJSPIsI.O 

6 0.377,  0,3FO,  0.3US,  0,333,  0.3??, 

7 0,310,  0,300,  0,?91,  0,?B5,  0,?83, 

TWTsO.O  ATMS?I=?,C 

« O.BOr,  O.iias,  0,U17,  r.39B,  r,-^73, 

9 0.3SO,  0,33?,  ^.3?1,  0.315,  0.311, 

TRTrO.o  AT>'’S?ls3,r. 

A o.soi,  O.aBP,  0./J93,  r,r4?i,  O,a03, 

B 0.377,  0.3^0,  0,351,  0.3t!3,  0,339, 

TPT  = ^ , 0 A7  ‘iSPiru  , 0 

C 0,50?,  0,075,  O.«ft0,  0.UU5,  0,U3?, 

n 0.910,  0,3B3,  0.350,  0.;<75,  0,355/ 

TBTsO.pS  AlPSaisO.O 

DATA  VAr75?/o.ia5,  o.i73,  O.io?,  0.?1?, 

1 0.??3,  0,29r,  0,?95,  0.?B0,  0,?50, 

2 0.?50, 

T9IsO.?5  AT^5?I=1,0 

3 0,?50,  0.?5?,  0,?55,  0.?5fe,  0,?57, 

9 0,?fe5,  0,259,  0,?70,  ri,27(i.  0,?55, 

TRT=o,?5  ATMS?I=2.0 

5 0.383,  o,3c;/j,  0,333,  0,317,  0,305, 

5 n.?Ql,  0.?B5,  0.?5P,  0.?51,  0.?51, 

TWT  = 0,?P  AT‘''S?I  = 3,0 

0.3B9,  0,371,  0.359,  0.390,  0,330, 

0.319,  0.303,  0,?90,  0.2B5,  o,?99, 

TPTsr.?5  Ar'''S?T  = 9,0 
9 0.909,  o,3Pg,  0.353,  0,375,  0,3P0, 

A 0,399,  0.^?5,  0.315,  0.305,  o.Tft/ 


i 


TABLE  XXVIII  (Continued) 


TRTsO  .*50  iT^'S?l=r,  .0 


1 

data  yAr783/0 . 070 , 0,138,  0.185, 

0.228, 

0,213, 

2 

0.235, 

0.238,  o.?u8,  o.2«7, 

TRIso.SC  4TMS2Iz1,o 

0.248, 

3 

0.179, 

0.203,  0.218,  0.229, 

0.238, 

4 

0.245, 

0.245,  0.245,  0.245, 

THI=0,50  ATMS2Is2.0 

0.245, 

5 

0,285, 

0.280,  0,248,  0,241, 

0.843, 

6 

0.248, 

0.252,  0.252,  0.250, 

TRTsO.50  AT‘:,S2I=3,0 

0.250, 

7 

0,305, 

0.280,  0,269,  0,280, 

0.257, 

8 

0.255, 

0,285,  0,255,  0.255, 

TRIeo,50  at'''S2I  = 4,0 

0.258, 

9 

0.319, 

0.300 , 0.288,  0,280, 

0.275, 

A 

O.270, 

O.270.  0.270,  0,270, 

TPTrl  ,00  ATM.S2Is0.0 

0,270/ 

1 

DATA  yAr7P-:4/0.000  , 0.1  1 7,  0.1  76, 

0.225, 

0.207, 

2 

0.239, 

0.243,  0.248,  0.246, 

TRjn.OO  AT^82Ts1,0 

0.248, 

3 

0.088, 

0.093,  0,108,  0.135, 

0.184, 

u 

0.208, 

0.222.  0.225,  0.230, 

TCT  = 1 ,00  AT^.S2Ts?.o 

0.231  , 

5 

0.168, 

0.187,  0.188,  0,189, 

0.170, 

b 

0,180, 

0,198,  0.208,  0,211, 

TBTsl.oo  AT  ►'.5  21  = 3.0 

0.214, 

7 

0.19ii, 

0.193,  o.ino,  0,190, 

O.190, 

8 

0,181, 

0.101,  0.198,  0.201, 

T9T  = 1.0  0 AT^'52T  = a,r 

0.203, 

9 

0,201  , 

0.200,  0.200,  0.197, 

0.104, 

A 

0.190, 

0.101,  0.  1 Of,  , 0.198, 

0,198/ 

TABLE  XXVIII  (Continued) 


CF^TfO  PF  PF  ltftIpG  S'IRPArFfSI.PF«sr>.lr) 


TPTsTiPFP  PATTP 

riATA  yAr7«)i/f',oro,  o.ps'L 

AT^>s?T  = A.SPFrr  RATir, 


0 . sn  0 , 


1 . nno , 


Tl'^FS  TA^-IGF^T  OF  MyrcWCPT  S^EFF 


1 

o.oor , 

1,000,  ? 

, 000  , 3 . 0 0 0 , 

9 , OOO , 

BART=ASPErT 

BATTfl  TTFFS  FFTA 

2 

0,000, 

0.500.  1 

. 000  , 1 .50'', 

?.ooo. 

3 

3.000, 

iJ  , 0 0 0 , 8 

,00'',  8.000, 

7 . 0 0 0 , 

y^'A^  = 8(‘PFACF 

rp  IN  pfpofI'T 

FEAN  GFOME'^RTC 

TPTsO.O 

ATNPPIzO.O 

a 

0.?F<0, 

0.2PP,  0 

.398,  0,900, 

0.950, 

S 

0.P8F, 

0.975,  0 

.980,  0.98?, 

0,985, 

TPIrO.O 

ATN.S?I  = 1 ,0 

b 

0.377, 

O.apo,  0 

.970,  0.970, 

0.970, 

7 

0.a70, 

0.981  , 0 

.985,  0.900, 

0.9P1  , 

TBTso.o 

atms?I=2.0 

P 

o.SOO, 

0,500,  0 

,500,  0,500, 

0.500, 

9 

O.FOO, 

0,800,  0 

.500,  0.500, 

0.500, 

TRTsO.O 

ATS)S?I  = 3.0 

A 

0.501 , 

0.518,  0 

.538,  0.536, 

0.538, 

« 

0.838, 

0.538,  0 

.538,  0.5P7, 

0,520, 

T P T * 0 . 0 

ATM5?I=9,0 

c 

0.50?, 

0,518,  0 

.839,  0.581, 

0.880, 

n 

0.580, 

0,580,  0 

.880,  0.580, 

0,858/ 

TPT=o,?5 

ATMS?!=O,0 

data  yAr7F?/o.]P5,  o 

.??5,  0.300, 

0.370  , 

1 

0 ,U50  , 

? 

'^.^J75, 

0.983,  0 

, 9 P 0 , 0 , 9 p 3 , 

0.9P8, 

TRTsO.PS 

AT‘i5?I  = l .0 

3 

0,?80, 

0.307,  0 

.390,  0.3P?, 

''.938, 

it 

0.P70, 

0,980,  0 

,988,  0,9pn, 

0 , 9P0 , 

TBIsO.?5 

AT^8?I=^.0 

S 

0.383, 

0,917,  0 

,995,  0.995, 

0,995, 

P 

0 .iia5. 

0 , upn , 0 

.500,  0.500, 

0,500  , 

tbIso,?5 

aTK'S?1  = 3.0 

7 

O.S'Jii, 

0,917,  0 

.991 , 0.985, 

0.980, 

8 

0 ,PPO  , 

0 .980  , 0 

.53?,  0.5?P, 

0.8??, 

TPIsO .?5 

A Tn?? 1 = 9 , 0 

P 

0 .AlOtJ  , 

0.917,  0 

.93?,  0.950, 

0,96P, 

A 

0.5?F, 

o.5?5,  0 

.8?8,  0.570, 

0.858/ 
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TPTro  ,8n 

A TMS?I  = 0 . 0 

hata  )(AC703/n,ri7n , 

0.18t,  0.257, 

0.358, 

1 

n.4J5,r , 

? 

0.4J7P, 

0.4188, 

0.4100,  C,aQi>, 

O.ap?, 

TWTrfl.50 

ATVS?I  = 1 .0 

3 

f'.17o, 

0,288, 

0.52P,  8,3fic, 

0 . '1 1 5 , 

a 

n.ijso, 

0,«75, 

0.880,  C.u«5, 

0 , a85 , 

TRT=0 ,S0 

ATMS?T=?.0 

5 

0.285, 

0.325, 

0,388,  0.^88, 

0.385, 

6 

0.^12?, 

0.4178, 

0 . 41 8 s , r , a p 0 , 

0 . apo , 

TRTsO ,80 

A T ^ s 2 I = 3 . 0 

7 

0.305, 

0.330, 

0.353,  0.371, 

0 . a 0 0 , 

8 

0.4138, 

0 . a8? , 

0.8)0,  0.508, 

0.808, 

T R T = n . 8 0 

4i  T582I=a.O 

Q 

0.31P, 

0.3?o, 

o.-^23,  0.331  , 

0.333, 

A 

.a6o , 

o.aoo , 

0.8)0,  0.5ao, 

0.538/ 

T P T = 1 . 0 0 

ATksPIsO.O 

DATA  VAr7P4J/O.OnO, 

0.000  , 0 . 330  , 

0 . a 1 2 , 

1 

0.41UT, 

2 

n,4ih8. 

o.aTi , 

0.4183  , 0,4185  , 

0 . 4188  , 

tptsi , on 

AT^'S2I  = 1 .0 

3 

0.088, 

o.?4'2. 

0.2P8,  0.388, 

0 , a 41  a , 

41 

n , a4ir> , 

0.881,, 

o.a88,  0.a70, 

0 . a70 , 

T R T r 1 . 0 

ATI  821  = 2,0 

5 

1 88, 

0.270, 

o.5ao.  0.37a, 

0,3PP, 

b 

0 . 4130  , 

0 . aaP , 

0 , ' 8 7 , C , a 8 8 , 

0 . 41  7 3 , 

TRT=1 ,00 

aT'',':?Is3.i, 

7 

0 . 1 04i  . 

''  .278, 

o.tap,  0.300. 

0 . 41  ) 2 , 

8 

n ,4i4i8  , 

0,858, 

o.a88,  o,a75. 

o.."78. 

T P T r 1 .80 

« T " 8 2 T = a . fi 

P 

0.201 , 

0.205, 

0,378  , (.,ii3a. 

0 . a 88  , 

A 

0 .4J7T, 

0 , a p ii , 

''.aPS, 

0 , aPa  / 
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6.  INLETS 

Three  types  of  inlets  are  used  for  the  airbreathing  configurations 
synthesized  in  the  CM-CGSM;  two-dimensional  side  mounted  (2  or  4)  inlets, 
and  a single  belly  line  mounted  inlet.  The  side  mounted  inlets  are  flat  sided 
and  mounted  symmetrically  with  the  internal  compression  surfaces  at 
the  top.  Belly  mounted  inlets  are  mounted  beneath  the  missile  body  with  the 
internal  compression  ramp  at  the  top.  A wedge  type  boundary  layer  diverter 
system  is  provided  for  both  inlet  types.  Figure  5 illustrates  the  inlet 
geometry  for  both  side  mounted  and  belly  mounted  inlets.  The  model  defines 
the  incremental  drag,  lift,  and  pitching  moment  associated  with  each  inlet 
type.  The  model  accounts  for  internal  (momentum)  lift,  but  does  not  account 
for  internal  drag  which  is  included  in  the  net  thrust  computation. 

6.  1 Drag 

The  incremental  drag  of  the  inlet  consists  of  a friction  drag, 
pressure  drag,  and  drag  due  to  lift. 

CDj  = CDj,  + CDp  + Cd. 

6.  1.  1 Zero  Lift  Drag 

Subsonic 

Zero  lift  drag  coefficient  of  the  inlet  subsystem  at  subsonic 
Mach  numbers  consists  o±  friction  drag  only.  The  net  friction  drag  incre- 
ment is  based  on  the  inlet  surface  area  minus  the  skin  area  of  the  main 
body  masked  by  the  inlet  and  fairing.  The  number  of  inlets  is  accounted  for 
by  the  factor  Njjsjl- 


^WET 

^REF 


N 


INL 


I 


The  skin  friction  coefficient  is  based  on  the  Prandtl-Schlicting  re- 
lationship for  tubulent  boundary  layers  corrected  for  compressibility  effects. 
Reynolds  number  is  based  on  the  total  inlet  length. 
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Transonic,  M = 1 

At  transonic  speeds  the  drag  consists  of  friction  drag,  forward 
fairing  drag,  and  aft  fairing  drag. 
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The  forward  fairing  drag  coefficient  is  based  on  a conical  fair- 
ing with  an  equivalent  conical  nose  fineness  ratio  of  3.  If  the  inlet  fairing 
extends  forward  of  the  non-tangency  point,  the  fairing  length  is  limited  to 
the  length  available  between  the  nose  tangency  point  and  the  inlet  leading 
edge.  The  diameter  of  the  equivalent  body  of  revolution  is  based  on  the  total 
projected  frontal  area  of  the  inlet.  Equivalent  cone  half  angle  is  based  on 
the  relation 


e = tan'^  (1.  /(FR  + FR) 
c 

The  forward  fairing  pressure  drag  at  M = 1 is  based  on  curve  fit  data  of 
Reference  3 presented  in  tabular  form  in  Table  VII. 

Afterbody  pressure  drag  is  based  on  an  equivalent  body  of 
revolution  defined  by  the  length  of  the  aft  fairing  and  the  equivalent  diameter 
discussed  above.  Empirical  data  of  Reference  4 is  used  to  estimate  the 
fairing  pressure  drag  at  M=1 

Cdpa  " 0.233/(L/D)2  . 0^^ 


Supersonic 

At  supersonic  speeds  below  the  ramjet  take-over  Mach  number 
the  inlets  are  covered  by  the  conical  fairings  discussed  above.  Zero-lift 
drag  coefficient  consists  of  the  skin  friction  drag  of  the  inlet,  pressure 
drag  of  the  forward  fairing,  and  pressure  drag  of  the  aft  fairing. 


5 
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Skin  friction  drag  is  computed  in  the  manner  described  in  the 
previous  sections.  Pressure  drag  on  the  forward  fairing  is  based  on  the  data 
of  Table  XII  for  an  equivalent  conical  body  of  revolution  with  a maximum  cross 
sectional  area  defined  by  the  total  projected  frontal  area  of  the  inlet.  After- 
body drag  is  computed  from  data  of  Table  XV  for  an  equivalent  body  of  revol- 
tion  defined  by  the  length  of  the  aft  fairing  and  the  equivalent  diameter  of  the 
inlet. 

At  supersonic  speeds  greater  than  the  ramjet  take-over  Mach 
number,  the  inlet  fairing  is  jettisoned  and  the  inlet  cowl  and  boundary  layer 
diameter  are  exposed.  The  boundary  layer  diverter  height  is  equal  to  75 
percent  of  the  turbulent  boundary  layer  thickness  at  the  inlet  leading  edge  at 
the  design  Mach  number  and  altitude.  A 10  degree  half  angle  wedge  is 
assumed  for  the  boundary  layer  diverter  system.  Diverter  system  projected 
area  consists  of  the  frontal  area  of  the  wedge  plus  an  additional  projected 
area  due  to  the  curvature  of  the  body.  Boundary  layer  diverter  wedge  drag 
is  based  on  data  of  Reference  2. 


The  term  (yTc  ) obtained  from  a curve  fit  to  data  of 

Figure  8,  page  17-9  of  Reference  2.  Mach  number  used  in  the  equation  is 
reduced  to  88  percent  of  the  free  stream  value  to  account  for  losses  in  the 
boundary  layer.  Friction  drag  of  the  diverter  system  is  based  on  the  wetted 
area  of  the  wedge  side  panels  and  the  friction  coefficient  is  based  on  Reynolds 
number  computed  at  the  diverter  station  based  on  body  length  to  that  point. 

Cowl  pressure  drag  is  defined  by  the  average  pressure  coefficient 
acting  on  the  cowl  projected  frontal  area.  Cowl  lip  pressure  coefficients  are 
computed  as  a function  of  Mach  number,  and  the  angle  of  the  oblique  shock 
waves  generated  by  the  e.xternal  compression  ramp  of  the  inlet  system. 
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Pressure  coefficient  on  the  cowl  lip  when  the  oblique  shock 
lies  inside  the  lip  is  based  on  the  following  relation: 

5 - 1 

^^LIP  ~ 3 ^.2 

If  the  shock  pattern  formed  by  the  multiple  compression  ramps  lie 
outside  the  cowl  lip,  a maximum  of  4 shocks  may  be  encountered  by  the  flow 
before  reaching  the  lip.  Oblique  shock  relations  derived  in  Reference  15 
are  used  to  define  flow  angle,  Mach  number  and  pressure  downstream  of 
each  shock  wave.  For  cases  when  an  oblique  shock  is  not  possible,  normal 
shock  relationships  are  used  for  the  case  of  an  oblique  shock 


Pressure  coefficient  is  computed  from  the  relation 

C = — - 1.  )/.  7 M^ 

P Pi 

Mach  number  downstream  of  an  oblique  shock  is  based  on  the 
shock  wave  angle  and  the  upstream  Mach  number 

-/)  ^ ^ ) 

{ /A/,*-  ^ 

If  the  flow  conditions  create  a normal  shock  pressure  ratio  across  the 

shock  wave  is  computed  by  the  following 


Pressure  coefficient  on  the  lip  is  computed  as  above.  The 
average  pressure  coefficient  acting  on  the  cowl  surface  is  assumed  to  be 
equal  to  1/2  of  the  pressure  coefficient  acting  on  a wedge  with  a constant  angle 
equal  to  the  lip  initial  angle.  Pressure  drag  coefficient  on  the  cowl  lip  is 
based  on  the  relation 

^°COWL  " ^^LIP  ^PROJ 
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Pressure  drag  coefficient  on  the  aft  fairing  of  the  inlet  is  based 


on  an  equivalent  body  of  revolution  defined  by  the  length  of  the  fairing  and 
its  equivalent  diameter.  The  data  of  Table  XV  is  used  for  fairings  with  an 
equivalent  fineness  ratio  of  at  least  2 to  1 . Inlet  fairing  pressure  drag  for 
inlet  fairings  with  fineness  ratio  less  than  2 to  1 is  computed  from  the  data 
for  conical  boattails  from  Reference  3.  If  the  pressure  drag  coefficient 
computed  for  the  short  afterbody  fairing  is  greater  than  the  full  base  drag 
on  the  inlet  projected  area,  base  drag  is  used  for  the  drag  of  the  afterbody 
fairing . 

6.1.2  Induced  Drag 

Induced  drag  is  equivalent  for  the  complete  configuration  as 
described  in  Section  5.1.2. 

6.  2 Lift  Curve  Slope 

The  two-dimensional  side  mounted  inlets  are  treated  as  a low 

aspect  ratio  lifting  surface  and  the  lift  curve  slope  is  computed  in  the  manner 

described  in  Section  5.  2.  Span  of  the  equivalent  lifting  surface  is 

equal  to  the  total  width  of  the  side  mounted  inlets.  Effective  area  of  the 

equivalent  lifting  surface  is  based  on  the  total  platform  area  of  the  inlet 

plus  one-half  of  the  platform  area  of  the  aft  fairing.  Aspect  ratio  of  the 

equivalent  surface  is  based  on  the  exposed  span  and  the  effective  area  of 

the  inlet.  A momentum  lift  term  based  on  slender  body  theory  is  used 

when  the  inlet  fairing  is  removed 

2. 0 A 

The  momentum  lift  term  is  the  only  lift  contribution  for  single 
belly  mounted  inlets. 

6.  3 Center  of  Pressure 

Center  of  pressure  for  the  side  mounted  inlets  is  located  at 
the  centroid  of  the  effective  planform  area  at  all  Mach  numbers.  Center 
of  pressure  for  the  belly  mounted  inlet  is  assumed  to  be  at  the  inlet  leading 
edge . 
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6.  4 


List  of  Symbols 

Symbol 

A 

c 

‘^Ddiv 

CD^ 

CDi 

CDp 


Cp. 


Description 
Inlet  capture  area 

Boundary  layer  diverter  pressure  drag 
coefficient 

Friction  drag  coefficient 
Induced  drag  coefficient 
Pressure  drag  coefficient 
Skin  friction  coefficient 

Compressibility  correction  factor 


FR 

L/D 

M 


Pressure  coefficient 

Lift  curve  slope  due  to  turning  of  the  flow 
at  the  inlet 

Fineness  ratio  of  inlet  fairing 

Length  to  diameter  ratio  of  inlet  aft  fairing 

Mach  number 


P2 

PI 

®REF 

^WET 

t/c 

^BT 


pressure  ratio  across  shock 
Reference  area 

Wetted  area  of  inlet  and  fairings 

Thickness  ratio  of  wedge  diverter 
Boattail  fairing  angle 
Oblique  shock  wave  angle 
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7.  EXTERNAL  BOOSTER 

The  external  booster  model  used  for  surface  launched  missiles 
is  shown  in  Eigure  8.  The  booster  case  is  a non-boattailed,  right  circular 
cylinder  with  a 25  degree  half  angle  nose  cone  with  a bluntness  ratio  of  25 
percent.  Booster  diameter  is  an  input  variable,  booster  length  is  a function 
of  the  propellant  required  to  achieve  a desired  flight  condition. 

7. 1 Drag 

Drag  coefficients  for  the  external  boosters  are  based  on  isolated 
body  coefficients  corrected  to  the  coefficient  reference  area  for  the  complete 
configuration. 

7.1.1  7ero-Lift  Drag 

External  booster  zero-lift  drag  consists  of  nose  cone  pressure 
drag,  base  drag,  and  skin  friction  drag.  Skin  friction  drag  coefficient  is 
computed  in  the  manner  previously  described.  Reynold's  number  computa- 
tions are  based  on  the  booster  length.  Nose  cone  pressure  drag  coefficient 
for  a 25  degree  half  angle  blunted  cone  is  presented  in  Figure  9.  Base 
pressure  drag  coefficient  is  based  on  the  base  pressure  coefficients  pre- 
sented in  Table  V.  Effective  base  area  is  50  percent  of  maximum  cross 
sectional  area  of  the  booster. 

7.1.2  Induced  Drag 

Induced  drag  is  evaluated  for  the  total  configuration  as  described 
in  Section  4.  1 . 2. 

7.  2 Lift  Curve  Slope 

Lift  curve  slope  for  the  external  booster  is  based  on  a correla- 
tion of  empirical  data  on  cone- cylinder  bodies.  Lift  curve  slope  for  the 
cone-cylinder  configuration  used  in  the  external  booster  model  is  shown  in 
Figure  10.  Lift  data  for  the  external  boosters  are  based  on  isolated  body 
coefficients  corrected  to  the  coefficient  reference  area  for  the  complete 
configuration. 

7.  3 Center  of  Pressure 

Center  of  pressure  data  for  the  external  boosters  is  based  on  a 
correlation  of  empirical  data  on  cone-cylinder  bodies.  Center  of  pressure 
of  the  isolated  booster  is  presented  in  Figure  1 1 as  a function  of  Mach  n\imber. 
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Lift  Curve  Scope 
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FIGURE  10  LIFT  CURVE  SLOPE 
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FIGURE  11  CENTER  OF  PRESSURE 
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ABSTRACT 


This  appendix  describes  the  SEATIDE  VEPiPER  sub- 
model. In  its  present  state,  VEHPER  computes  the  trajectory  of  a 
cruise  missile  powered  by  rocket , tur  bojet  or  ramjet  propulsion  systems. 
The  description  includes  math  models  for  the  vehicle  characteristics  and 
its  geophysical  environment,  numerical  techniques,  and  a top-level 
block  diagram. 
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1.  INTRODUCTION 


The  role  of  the  VEHPER  submodel  is  to  compute  the 
performance  of  atmospheric  cruise  missiles  synthesized  by  the  SEATIDE 
Concept  Generation  and  Screening  Model.  The  submodel  is  capable  of 
analyzing  a wide  variety  of  missions  and  configurations.  The  vehicle  may 
be  surface  or  air-launched,  and  it  may  fly  high  altitude,  low  altitude,  or 
split-level  missions.  Propulsion  systems  represented  include  a rocket 
booster  and  throttleable  rocket,  ramjet,  and  turbojet  sustainers. 

The  basic  approach  to  computing  the  performance  of  the 
missile  is  to  numerically  integrate  two  degree-of-freedom,  point  mass 
equations  of  motion.  The  integration  subroutine  utilizes  a fourth  order, 
variable  step  size  Runge-Kutta  method  which  controls  the  integration 
step  size  to  maintain  an  approximately  constant  truncation  error.  Seven 
first-order  differential  equations  of  motion  are  integrated  simultaneously. 
Variables  integrated  include  flight  path  variables,  weight,  and  ideal 
velocity. 

The  trajectory  is  divided  into  phases  for  computational  and 
user  convenience.  At  the  start  of  each  phase,  the  vehicle  characteristics 
and  its  flight  path  control  method  may  be  specified  in  order  that  the 
vehicle  can  fly  through  a sequence  of  configurations  and  maneuvers. 

Figure  1 illustrates  a typical  sequence.  Initial  conditions  define  a low- 
j altitude  air  launch.  In  the  first  phase,  the  missile  is  accelerated  to  a 

j desired  Mach  number  by  a rocket  booster.  In  the  second  phase, 

I aerodynamics  are  re-defined,  weights  associated  with  the  boost  phase 

' are  dropped,  and  the  sustainer  propels  the  missile  during  a 

climb.  The  third  phase  is  a maneuver  for  leveling  the  climb  path  at  the 
cruise  altitude.  Phase  4 is  the  cruise  phase.  Phase  5 is  a descent. 

Phase  6 is  a pull-out  at  a low  altitude.  The  7th  and  final  phase  is  a low- 
altitude  run-in  to  the  target. 
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From  a propulsion  standpoint,  a phase  may  be  one  of 


throe  types.  Tlie  phase  may  be  unpovvered,  (thrust  and  fuel  flow  are  zero) 
and  the  vehi(.-le  is  in  "coasting"  flight  affected  only  by  aerodynamic  and 
gravity  forces.  A coast  phase-  might  be  utilized  to  provide  a realistic  time 
delay  l:>etween  Ijooster  burnout  and  air-breathing  sustainer  ignition  or  an 
unpowered  clescent  to  target  encounter.  A second  type  of  propulsion  in  a 
phase  is  boost.  A boost  phase  is  one  in  which  an  unthrottled  rocket 
propels  the  v'ehicle.  A boost  phase  would  typically  be  used  to  accelerate 
a ramjet  from  a subsonic  (or  zero)  velocity  to  a supersonic  ramjet  take- 
over Mach  number.  The  third  category  is  one  in  which  the  vehicle  is 
propelled  by  a sustainer.  The  sustainer  is  characterized  as  being  capable 
of  throttled  operation  in  order  to  satisfy  a desired  velocity  requirement, 
e.  g.  , a cruise  Mach  number.  At  the  present  state  of  development,  the 
sustainer  ma>'  be  either  a rocket,  ramjet,  or  turbojet.  For  the 
duration  of  a phase,  the  sustainer  may  be  operated  i.i  one  of  three  modes; 
throttleable,  maximum  thrust,  or  minimum  thrust. 

Some  important  optional  capabilities  have  been  added  to 
the  basic  trajectory  program  to  relieve  the  user  of  tedious  "hunt  and  try" 
procedures  for  finding  acceptable  flight  paths. 

In  assessing  a missile  design,  the  analyst  must  be  able 
to  compute  its  maximum  range.  For  the  mission  profile  of  Figure  1, 
maximum  range  can  be  achieved  by  increasing  the  length  of  the  cruise 
leg  of  the  mission  until  luel  is  exhausted  at  the  end  of  the  low  altitude 
run-in.  The  V'EHPLR  submodel  allows  the  analyst  to  designate  a 
particular  phase  as  a variable  length  cruise  phase,  and  the  computer 
program  will  perform  an  iterative  search  for  the  cruise  leg  length  which 
exhausts  fuel  at  the  end  of  the  last  phase. 

Another  requirement  placed  on  the  routine  is  that  the 
climb  (o  altitude  be  performed  efficiently,  as  experience  has  shown  that 
a climb  path  angle  which  is  too  great  or  too  small  can  hav'e  a drastic 
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effect  on  total  range.  Rather  than  require  the  user  to  guess  at  the  proper 
climb  angle,  provision  has  been  rnacle  to  generate  a climb  schedule,  con- 
sisting of  a Mach  number  versus  altitude  table.  Just  before  the  climb 
phase  is  begun,  the  table  is  constructed  by  searching  for  the  best  climb 
Speed  for  each  of  several  different  altitudes,  erod  ynam  ics  and  propul- 
sion characteristics  are  utilized  to  determine  excess  thrust  at  each 
velocity  and  altitude  combination  considered. 

A third  capability  designed  to  ease  the  user's  task  is 
the  incorporation  of  a novel  level-off  navigation  law.  During  climbs  or 
descents  to  a cruise  altitude,  the  level-off  maneuver  must  be  initiated  at 
a point  on  the  path  which  allows  for  the  time  to  change  the  path  angle  to 
zero  while  the  cruise  altitude  is  being  achieved.  The  altitude  change 
during  the  level-off  maneuver  can  be  several  miles  for  a Mach  4 ramjet. 
The  point  of  level-off  initiation  could  be  found  through  the  "hunt  and  try" 
process  of  running  several  level-off  trajectories  starting  at  different 
times  on  the  climb  or  descent  path.  The  VEHPER  eliminates  the  need 
for  such  iterative  techniques  by  employing  a simple,  approximate 
analytical  expression  relating  level-off  load  factor,  altitude  error  with 
respect  to  desired  cruise  altitude,  and  vertical  speed.  The  climb  or 
descent  phase  is  stopped  when  the  load  factor  required  to  level  off  at  the 
cruise  altitude  becomes  equal  to  the  level-off  load  factor  desired  by  the 
user.  From  the  climb  or  descent  stopping  point  until  reaching  the 
cruise  altitude,  the  path  is  controlled  by  the  level-off  navigation  law. 

At  each  integration  step  along  the  trajectory,  the  load  factor  required  is 
computed  from  the  altitude  error  and  vertical  speed.  A constant 
velocity  during  the  level-off  will  result  in  a near-constant  load  factor. 


2.  FLIGHT  PATH  EQUATIONS 


EQUATIONS  OF  MOTION 


Two  dimensional  flight  path  equations  of  motion  define  the 
missile  s veloc  ity  rate  and  path  angle  rate  with  respect  to  a non- rotat ing , 
spherical  earth  model; 


V = Tco<^o<-»n)  -D  - ,3  *1^  J,*  (f+/*ec') 


(2-0 


y = ^ (raa/$ec‘)  C2-2) 


Rates  of  change  of  altitude  and  down  range  distance  are  respectively 

X ~ V fiw  y (ft/set.) 

R = V Rc  COi  t ( ft/sec) 

r 

Initial  conditions  for  the  differential  equations  of  motion  are: 


V(o)  = 

MAX  (l/o,  0.1) 

y(o)  = 

It,  * Tf/l80 

X(o)  = 

X, 

(fO 

R(0)  = 

R.  * 4074.115 

CfO 

where  Voi  ^ ‘iltf  R c are  inputs.  Path  variables  used  in  the 
equations  of  motion  are  illustrated  in  Figure  2 and  are  defined  as  follows: 
V = velocity  (ft/sec) 

T = net  thrust  force  (Ibf) 

0^  = angle  of  attack  (rad) 

D = aerodynamic  drag  force  (Ibf) 

KVI  = mass  (slugs) 

= gravity  force  per  unit  mass  (ft/sec^) 
if  = flight  path  elevation  angle  (rad) 

L = aerodynamic  lift  force  (Ibf) 

f = geocentric  radial  distance  to  vehicle  (ft) 

H = thrust  cant  angle  for  boost  phases  (radians) 

Vehicle  mass  is  determined  by  integrating  its  time  rate, 

M*"  = - AO'f  C2-s) 


and  juiXf  is  the  fuel  flow  rate  computed  by  the  propulsion  math  model. 

The  missile  angle  of  attack  is  the  control  variable  used  to 
determine  the  flight  path;  its  value  is  constrained  to  an  input  maximum 
value  such  that 


kl  ± » Tr/l80 

where  input  maximum  in  degrees.  Subsequent  paragraphs  will 

indicate  the  several  options  for  computing  angle  of  attack. 

The  missile  flight  path  is  computed  by  simultaneously 
integrating  five  first-order  differential  equations,  Eq.  (2-1),  (2-2),  (2-3), 
(2-4),  and  (2-5).  The  numerical  integration  procedure  is  a fourth  order 
variable  step  size  technique  described  in  Section  3.  0, 
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FIGURE  2 FLIGHT  PATH  VARIABLES 


PATH  CONTROL  OPTIONS 


2.  2 

The  missile  angle  of  attack  is  the  final  output  of  each  of 
the  flight  path  control  options.  The  cruise  and  level-off  options  also 
determine  thrust  required. 

2.2.1  Option  _Li  Jf(Jt)  Path  Angle  Control 

An  input  table  of  flight  path  angle  versus  time  determines 
the  desired  angle  of  attack.  The  commanded  path  angle  is 

a rU) 

where  is  bounded  by  tabular  values,  ^(^4 1 and  , and  iL  is 

likewise  bounded  by  tabular  values  and  • .A  computed  path 

angle  rate  is 

The  angle  of  attack  required  to  achieve  F is  determined  by  an  iterative 
solution  of  the  If  equation  of  motion,  Eq.  (2-2).  Subroutine  XALPHl 
performs  this  function.  The  variation  in  thrust  with  angle  of  attack,  if  any, 
is  considered  in  XALPHl, 

2.2.2  Option  2:  Not  Used 

2.2.3  Option  3;  9(t)  Attitude  Control  vs . Time 

Angle  of  attack  is  computed  directly  from 

o(  s fi(t)  - t (t^A) 

2.  2.  4 Option  ± oU)  Attitude  Control  vs.  Altitude 

Angle  of  attack  is  computed  directly  from 

(faa) 

2.2,5  Option  5;  Attitude  Control  vs . Mach  No. 

Angle  of  attack  is  computed  directly  from 

oC  = © (m)  - (ra4) 
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2.2.6 


2.  2.  7 


Option  6:  CX,Ct)  Anj>le  of  Attack  Control  vs.  Time 
Angle  of  attack  is  determined  by  a table  look-up 


ei  - oi 


(t)  (rad) 


Option  7:  o((A)  Angle  of  Attack  Control  vs.  Altitude 

Angle  of  attack  is  determined  by  a table  look-up 

oi-et^CJc)  (rad) 


2.2.8  Option  8 : Angle  of  Attack  Control  vs.  Mach  No 

Angle  of  attack  is  determined  by  a table  look-up 

= o<  (k)  (rad  ) 

2.2.9  Option  9:  Nj(it)  Factor  Control  vs.  Time 

Load  factor  is  determined  by  a table  look-up 

= (g*0 

Path  angle  rate  as  a function  of  load  factor  is 


V (raetisat) 

The  angle  of  attack  to  achieve  the  path  angle  rate  is  computed  in 
XALPHl  which  iteratively  solves  Eq.  2-2  for  o(  . 

2.2.  10  Option  10:  Load  Factor  Control  vs.  Altitude 

Load  factor  is  determined  by  a table  look-up 

= Cg'«) 

Path  angle  rate  and  angle  of  attack  are  computed  as  in  Option  9. 
2.2.  11  Option  1 1 : N Load  Factor  Control  vs.  Mach  No . 

Load  factor  is  determined  by  a table  look-up 

= (&'») 

Path  angle  rate  and  angle  of  attack  are  computed  as  in  Option  9. 


^ = [ COS  j 


s ubroutine 

i 


A 
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2.2.  12 


Option  12:  Mach  Number  Climb  or  Descent  Schedule 

This  is  the  primary  path  control  option  for  performing  the 
climb  and  descent  phases  on  a mission  profile.  The  table  is 

optionally  input  to  VEHPER  or  internally  computed  within  the  VEHPER 
MAINS  subroutine.  The  engine  is  assumed  to  be  operating  at  maximum 
thrust;  hence  the  schedule  is  maintained  by  controlling  the  path  elevation 
angle. 

The  path  control  equation  to  achieve  the  desired 

schedule  is 

/=  Kv  + ky  ( IT'*- jr)  (raJi/itc') 

where 

are  input  gains 

is  the  desired  velocity  (ft/sec); 

V*5M(iL)*0t  (ft/sec) 

^ = sound  speed  at  altitude  (ft/sec) 

is  a computed  desired  path  angle  (rad) 

The  desired  path  angle  is  based  on  the  slope  of  the  Mai  curve,  thrust, 
and  drag; 


1^*  a sm  T*co>  (r*4) 

is  the  angle  of  attack  on  the  desired  path  (radians) 
is  the  drag  on  the  desired  path  (Ibf) 
is  the  thrust  on  the  desired  path  (Ibf) 

= rate  of  change  of  velocity  on  the  schedule 

Since  there  are  two  unknowns,  if*  and  O^*,  in  Eq.  (2-6)  an  additional 
consideration  must  be  introduced  to  permit  the  solution.  An  assumption 
that  IT  = 0 yields 

toi  If*/ (t*  Ci.^  a- S)  (rad)  (l-7) 


0' 

T‘ 

die 
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Equations  (2-6)  and  (2-7)  are  solved  simultaneously  for  IT  and  o(  in 

an  iterative  procedure. 

Z.2.13  Option  13:  Cruise 

During  crviise  flight  the  vehicle  is  flown  at  constant  altitude 

(if  the  angle  of  attack  constraint  permits),  and  sustainer  engine  thrust  is 
controlled  to  satisfy  the  desired  cruise  Mach  number. 

The  concept  of  accelerations  normal  and  tangent  to  the 
flight  path  are  introduced  for  the  cruise  path  control  option  (and  also  the 
Option  14  level-off  path  control).  Accelerations  normal  and  tangent  to  the 
path  are  respectively 


f T avAot  U 4.  ( 1/*--  < 

(i-a) 

1 VIA  r 

' \ 

P 

-1 

II 

!<• 

II 

T coioit  “0  - 

(4'0 

(z-3) 

vv\ 

1 

During  cruise  flight, 

an  =0 

(c»'j) 

(2-10') 

ax  . («‘a- v)  C4'd 

KtnO  Is  a gain  on  the  velocity  error  (G’s/ft/sec) 

W**  is  the  desired  (cruise)  Mach 
OL  is  the  speed  of  sound  (ft/sec) 

V is  the  missile  velocity  (ft/sec) 

Subroutine  XALPH2  solves  equations  (2-8)  and  (2-9)  for 
required  thrust  T and  angle  of  attack  o<  which  satisfy  cruise  conditions 
specified  by  equations  (2-10)  and  (2-11).  If  thrust  required  is  not  available 
for  the  sustainer  engine,  then  maximum  thrust  is  commanded,  and  angle  of 
attack  is  recomputed  based  on  maximum  thrust.  If  angle  of  attack  required 
for  Oln  = 0 exceeds  the  maximum,  the  vehicle  is  flown  at  maximum  angle  of 
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attack  and  a descent  to  a lower  altitude  will  occur.  Maximum  angle  of 
attack  flight  continues  until  the  flight  path  angle  becomes  zero,  and 
thereafter  the  vehicle  is  flown  at  O-M  = 0. 


1 

1 


2.2.  14  Option  14;  Level-Off 

The  level-off  maneuver  is  used  as  a transition  maneuver 
between  a climb  or  descent  phase  and  a constant  altitude,  (zero  path  angle) 
cruise  phase.  A load  factor  is  computed  which  will  cause  the  path  angle  to 
become  zero  at  the  same  time  that  cruise  altitude  is  attained.  The  sustainer 
propulsion  system  is  throttled  to  maintain  a desired  Mach  number. 

The  level-off  navigation  law  relates  acceleration  normal 
to  the  path,  QLn  , to  altitude  error  and  altitude  rate: 

= "A  I 4:.  1 

2^  I 

= altitude  error;  cruise  less  current  altitude  (ft) 

The  singularity  at  cruise  altitude  is  circumvented  by  holding  (!»•  constant 
during  the  last  20  feet  of  altitude  change.  The  law  yields  a constant  normal 
acceleration  if  the  missile  velocity,  V/  , is  constant. 

Tangential  acceleration  during  level-off  is  computed  from 
the  same  expression  used  for  cruise: 

Q-t  * (h^a.  - v)  C&'*) 

Just  as  in  the  cruise  option,  subroutine  XALPH2  solves 
for  angle  of  attack  and  thrust  which  satisfy  the  required  normal  and 
tangential  acceleration.  If  thrust  required  is  greater  than  maximum  avail- 
able, the  sustainer  is  operated  at  maximum  thrust. 

2,3  AERODYNAMICS  FORCES 

Aerodynamic  lift  and  drag  forces  on  the  vehicle  are 

respectively 

L = Ci.«^  ot  S ( lUf'i 

D=  C Cd,  ^ (Ibf) 
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L ■ ^ 


where 


2.4 


s 

Cdo 


is  the  slope  of  the  lift  curve  and  is  an  input  tabular 
function  of  Mach  number  (radians 
is  the  angle  of  attack  (radians) 
is  the  dynamic  pressure;  % ~ x ^ ^ ^ ’ 

^ is  the  atmospheric  density  in  slugs/ft 
(Ibf/ft^) 

2 

is  the  aerodynamic  reference  area  (ft  ) 
is  the  zero  angle  of  attack  drag  coefficient; 
an  input  tabular  function  of  Mach  number  and 
altitude. 


PROPULSION 


At  the  present  State  of  development,  the  VEHPER  submodel 
considers  rocket,  turbojet,  or  ramjet  propulsion.  The  rocket  may  be  either 
a constant  thrust  booster  or  a throttleable  sustainer.  The  ramjet  and  turbojet  sus 
tainers  may  be  operated  in  a throttleable  mode  or  at  maximum  or  minimum  thrust 
2.4.1  Rocket  Booster 

Both  liquid  and  solid  rocket  thrust  and  fuel  flow  are 
modeled  with  the  following  relationships: 


Mt  = T/lt» 


( IWm  / s*c) 


T 

At 

I»r 


is  net  thrust  ( Ibf ) 

is  the  input  constant  vacuum  thrust  (Ibf) 

2 

is  the  rocket  nozzle  exit  area  (ft  ) 

2 

is  the  ambient  atmospheric  pressure  (Ibs/ft  ) 
is  fuel  flow  rate  (Ibm/sec) 

is  the  input  constant  vacuum  specific  impulse 
(Ibf  sec/lbm) 


2.4.  2 


Rocket  Sustainer 


During  throttled  operation,  thrust  required,  T*  , is 
computed  to  satisfy  trajectory  requirements  during  each  integration  step. 
The  vacuum  thrust  required  is 

W = T*  4 (Itf) 

and  is  bounded  according  to 

( t TyfKt.  * 

where 


^ (Tvkc)(^(^5(  are  input  maximum  and  minimum  values 
of  vacuum  thrust  (Ibf) 

The  throttle  fraction  is  the  ratio  of  the  thrust  increment  used  to  the  thrust 
increment  available: 


R = 


~ CTv^c)mim 
“CTvAt)HXH 


Fuel  flow  is  computed  from  a vacuum  specific  impulse  which  is  an  input 
tabular  function  of  throttle  fraction: 

TspCa') 

Net  thrust  force  is 

T = - Ae  p,  ( IWf) 


2.4.3  Ramjet  Sustainer 

The  source  of  ramjet  performance  data  is  subroutine 
PROPll.  The  equations  are  based  on  compressible  flow  laws  of  gas 
dynamics  and  thermodynamics.  The  math  models  employed  by  PROPll  are 
those  detailed  for  PROPl  in  Appendix  F.  Because  the  computations 
produce  lengthy  computer  execution  times,  the  representation  of  ramjet 
characteristics  has  been  linearized  to  reduce  the  number  of  times  the 
PROPll  subroutine  must  be  executed.  In  the  course  of  trajectory 
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calculations,  linear  relationships  for  ramjet  performance  characteristics 
are  utilized  in  the  equations  of  motion.  The  PROPllsubroutine  is  used 
only  to  produce  partial  derivatives  for  the  linear  relationships.  Partial 
derivatives  are  computed  by  finite  differences. 

Ramjet  performance  linearization  is  based  on  the  assumptions 
that  performance  characteristics  are  linear  functions  of  flight  conditions 
within  a flight  condition  region.  A region  is  established  at  the  beginning 
of  ramjet  powered  flight  and  thereaf.e  is  re-established  each  time  a 
flight  condition  exceeds  a bound  of  the  previous  region. 

Maximum  or  minimum  thrust  engine  operation  is  utilized 
respectively  for  climbs  or  descents.  Ramjet  thrust  coefficient,  specific 
fuel  flow,  and  combustion  temperature  are  represented  as  linear  functions 
of  angle  of  attack, altitude,  and  Mach  number  Thus  for  maximum  or 
minimum  thrust  flight  (designated  by  an  input  flag). 


T 


+ ii?  (rtk-ot.)  + SCp  Ci-X.) 


+ aCF  a S 

SM  J ® 

Mff  ~ I Sf«  + (X--R.e') 

^ a d 
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where  performance  characteristics  are 


T net  thrust  (Itjf) 

Cf  'let  thrust  coefficient 

fuel  flow  (Ibm/sec) 

Sp  specific  fuel  flow 

o 

Tt^  combustion  temperature  ( R) 

and  flight  conditions  are 

0<w  angle  of  attack  (rad) 

altitude  (ft) 

M Mach  number 

Nominal  values  of  the  performance  characteristics  and  flight  conditions 
(values  at  the  "center"  of  the  flight  condition  region)  are  denoted  by  the 
"0"  subscript. 

Throttled  engine  operation  is  utilized  during  cruise  or 
during  the  level-off  maneuvers  to  achieve  a desired  Mach  number.  For 
this  case,  the  required  thrust  coefficient  is  an  independent  variable,  and 
we  have 


f 


4 £Sf  Cto') 

acf  J 

+ ( a. - 

sir 

+ ( Cf  “ Cf,') 

SCf 

R 


r 


t 

I 


t 

I 


I 


i 
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where  Cp  is  the  desired  thrust  coefficient  as  determined  by  trajectory 
considerations,  and  Cp,  is  the  nominal  values  for  the  established  flight 
condition  region.  Bounds  on  desired  thrust  coefficient,  which  are  used  to 
represent  maximum  and  minimum  thrust  available,  are 

0 f Cp  i Cp^,^, 

The  maximum  thrust  coefficient  must  also  be  represented  as  a function  of 
flight  conditions; 


arA 

The  zero  lower  limit  on  net  thrust  coefficient  is  generally  not  representative 
of  the  minimum  available;  negative  values  are  usually  achiev:.ble  during 
ramjet  operation  at  minimum  fuel  flow.  However,  the  difference  in  the 
zero  limit  and  the  achievable  negative  value  should  not  affect  the  trajectory 
to  a significant  degree  for  foreseeable  applications.  The  use  of  a zero 
limit  serves  to  reduce  computer  execution  time  and  to  ease  PROPl  con- 
vergence difficulties  associated  with  negative  thrust  coefficients. 

2.4.4  Turbojet  Sustainer 

The  calculation  of  thrust  and  fuel  flow  for  a turbojet  sustainer 
entails  the  computation  of  pressure  recovery  and  corrected  airflow  which 
satisfy  both  the  inlet  and  the  engine  requirement.  For  a given  engine 
design,  a specified  turbine  inlet  temperature,  and  flight  conditions  (Mach 
number,  altitude,  and  angle  of  attack),  the  relationships  between  pressure 
recovery  , the  corrected  airflow  are  illustrated  in  the 

sketch. 
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It  is  seen  that  if  flight  conditions  and  the  turbine  inlet  temperature  result  in 
subcritical  corrected  airflow,  the  inlet  can  deliver  any  corrected  air  flow 
less  than  the  critical  value  at  the  design  pressure  recovery.  (Critical  flow 
conditions  result  in  a normal  shock  at  the  inlet  throat;  subcritical  flow  means 
that  the  shock  is  upstream  of  the  throat;  and  supercritical  flow  means  that 
the  shock  is  downstream  of  the  throat.  ) If  the  engine  demands  a super- 
critical corrected  airflow,  then  the  engine  and  inlet  relationships  must  be 
solved  simultaneously  for  pressure  recovery  and  corrected  air  flow.  The 
algorithm  for  the  engine  corrected  air  flow  requirement  as  a function  of 
pressure  recovery  is  embodied  in  GENENG,  a subroutine  version  of  the 
computer  program  documented  in  Reference  2.  The  inlet  algorithm  is  con- 
tained in  subroutine  INLET  and  is  described  in  Appendix  E.  Subroutine 
TJPER  controls  the  calculations  performed  by  INLET  and  GENENG. 

2.5  GRAVITY 

The  gravitational  force  per  unit  mass,  acting  at  the  vehicle 
mass  center  and  directed  toward  the  center  of  the  spherical  earth  model,  is 
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is  zero-altitude  standard  gravity; 

32.  1 7405  ft/sec^ 

is  the  radius  of  the  spherical  earth  model; 
20.  9 X 10^  ft 
is  vehicle  altitude  (ft) 

2.6  ATMOSPHERIC  PROPERTIES 

A time  and  location  invariant  U.S,  Standard  1962  model 
atmosphere  is  used  to  define  atmospheric  properties  including  density, 
pressure,  temperature  and  speed  of  sound.  The  standard  atmosphere 
equations  are  contained  in  Reference  (1). 

2.  7 ROCKET  ENGINE  IDEAL  VELOCITY 

Rocket  engine  ideal  velocity  is  computed  in  the  course  of 
trajectory  computations  in  order  that  velocity  losses  can  be  quantified. 
Ideal  velocity,  the  velocity  which  the  vehicle  would  attain  in  the  absence  of 
atmospheric,  gravity,  and  turning  losses,  is 


k 


where 


(ft/5«c) 


respectively  the  rocket  engine  ignition 
and  burnout  times  (sec) 
is  vacuum  thrust  (Ibf) 
m is  vehicle  mass  (slugs) 

Total  velocity  loss  may  be  computed  from 


where  V (itO  “ V 

engine  operation. 


is  the  actual  velocity  change  during  rocket 


B-21 


2.8 


CLIMB  SCHEDULE  GENERATION 


A climb  schedule  can  be  generated  just  before  the  start 
of  a climb  phase  if  the  user  selects  this  optional  capability.  The  climb 
schedule,  which  consists  of  a table  of  Mach  number  versus  altitude,  is  part 
of  the  phase  initialization  process  accomplished  in  the  MAINS  subroutine. 

The  principal  relationships  utilized  are  the  equation  of 

motion,  Eq.  (2-1)  and  (2-2).  The  assumptions  are  made  that  missile 

• • 

velocity  rate  V , path  angle  rate  , and  thrust  cant  angle  •(  are  zero 
and  that  mass  is  constant  during  the  climb.  The  forms  of  the  equations  of 
motions  utilized  are  obtained  by  solving  Eq.  (2-1)  for  sin  If  and  Eq,  (2-2) 
for  o4  : 

•m  f “ L T "D 

% 

at,  - W'  ( ( T + Cu^  ^ s')  (rtJ)  (a -13^ 

Maximum  rate  of  climb  is  obtained  by  maximizing  the 
quantity  V at  each  altitude  of  an  array  of  altitudes  encompassing  the 

expected  climb  path.  At  each  velocity-altitude  combination, Eqs.  (2-12) 
and  (2-13)  are  solved  simultaneously  for  the  angle  of  attack  and  path  angle 
by  means  of  an  iterative  successive  substitutions  technique.  The  search 
for  climb  speed,  , is  limited  according  to 

Vmiw  - i Vmkh 

where 

V/ftv,  = the  climb  speed  (ft/sec) 

sec) 

is  an  input  reference  Mach  number 
for  the  phase 

QL  is  the  local  speed  of  sound  (ft/sec) 
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the  altitude  are  computed  and  stored  in  the  climb  schedule  table. 

Because  of  the  assumptions  made  in  formulating  the  climb 


schedule  generation,  limitations  on  its  use  should  be  observed.  Obviously 
a vehicle  with  a high  thrust  to  weight  ratio  would  not  be  consistent  with 
the  assumption  that  V*0,  and  in  fact  thrust-to-weight  ratios  greater 
than  unity  can  cause  Eq,  (2-12)  to  yield  a value  for  sin  which  is  greater 
than  one,  a mathematical  impossibility.  The  assumption  of  constant  mass 
is  not  normally  significant,  but  it  could  be  significant  if  a large  percentage 
of  the  initial  mass  is  consumed  during  the  climb  (e.  g.  a low  thrust  to 
weight  ratio  rocket  with  an  extended  climb)  or  if  weight  is  jettisoned  during 
the  climb.  For  vehicles  where  these  modeling  limitations  preclude  internal 
climb  schedule  generation,  a climb  schedule  (Mach  No.  vs.  altitude)  can  be 
provided  by  the  user  through  Option  12  path  control  (see  Section  2.  2.  12). 
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3,  NUMERICAL  METHODS 


3.  1 RUNGE-KUTTA  INTEGRATION  WITH  ERROR  CONTROL 

The  Runge-Kutta  formula  used  is  of  fourth-order  accuracy 
in  step  size  h.  It  is  of  the  form 


where 


Xj  ^ - Xi  " X)  — g ( K|  + Z ♦ 2 Kj  4 ^4") 


a dependent  variable 

increment  in  the  dependent  variable 

increment  in  the  independent  variable  t 


K|  — Xi  C ^ X 1^ 

Kj  = AtX»(*,+  x,4«') 

Kj  * At  Xt  ( X(  + 

2 2 ' 

K4  S , Xi  + 

A lower-order  formula  may  be  found  by  utilizing  the  three 
derivatives  at  2 s , 3fc|  , and  X*  . If  and  Ats  ti-2, , 

the  following  Lagrangian  interpolation  formula  gives  the  derivative  at  any 
time  2 S 

X » X>  (t  _ X|  (t-t,)(2~ti) 

4 Xi  tt") 
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Integration  of  this  equation  from  to  yields 


4 ^ 2 -t.1  + _Al ^ X 


The  difference  in  t he  increments  over  the  interval  between  the  Runge- 

Kutta  scheme  and  the  low-order  scheme  may  be  divided  by  a nominal  value 
of  the  dependent  variable  to  obtain  the  relative  error  Thus, 


u = 


xt  - xr 


The  error  is  expected  to  vary  as  approximately  the  fifth 
power  of  , which  leads  to 

^ * A X* 

(where  A is  a suitable  coefficient)  or  in  the  logarithmic  form 

lo^  S ® A*  r 

where 

a'  * a 

Let  it  be  assumed  that  A'  will  vary  linearly  with  the  variable  of  integra- 
tion, Then  A'  at  a time  corresponding  to  can  be  found  from  A'  at  two 
previous  points  and  t,x 

A'a  = a'^  + - a!  iU-tC) 

and  if  s ~ it.  i 

Ai  = A»  + (Ax  - A>  ) 

and  on  this  basis  would  be  predicted  to  be  Aj  + S 
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It  is  desired  that  should  approximate  S , the  reference  error;  therefore. 

Each  dependent  variable  has  an  associated  relative  error  and  would  lead  to 
computation  of  a different  step  size  for  each  variable;  however,  the  maxi- 
mum relative  error  of  all  variables  may  be  selected  for  % . Obviously, 

inaccurate  predictions  of  step  size  can  occur  when  the  maximum  relative 
error  shifts  from  one  variable  to  another  or  when  any  sudden  change  occurs. 
When  a step  size  produces  an  excessively  large  error  Ci  ^ ^ 

reduced  step  size  must  be  used.  It  may  be  obtained  from  the  reference 
error  ^ as 

3.  1.  1 Starting  the  Integration 

The  Runge-Kutta  scheme  is  simple  to  start,  since  integration 
from  X*»  to  requires  no  knowledge  of  X prior  to  Xn.  Since  the  error 

control  coefficient  A has  no  value  at  0,  however,  a prediction  of  the  second 
step  size  is  difficult.  To  overcome  this  difficulty,  two  equal  size  first  steps 
may  be  made  before  checking  the  error.  The  A for  the  first  step  may  be 
arbitrarily  set  equal  to  the  A for  the  second  step  so  that  Jif  may  be  pre- 
dicted. The  low-order  integration  scheme  equation  in  this  case  becomes, 
with  JLt  , 

xT=  +4X,  + X,) 

3,1.2  F allures 

Should  two  consecutive  predictions  of  the  same  step  fail  to 
produce  an  error  S less  than  ^ return  to  the  starting  procedure 

will  be  made  with  a third  prediction  on  step  size,  which  is  no  larger  than 
one-half  of  the  second  estimate.  The  step-size  control  described  here 
will  operate  stably  with  nearly  constant  error  per  step  only  for  a well- 
behaved  function.  For  most  problems  it  will  repeat  a step  occasionally  to 
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reduce  a large  error,  and  on  sharp  corners  it  will  restart.  This  action 
is  not  regarded  as  objectionable.  The  objective  is  to  attain  a desired 
level  of  accuracy  with  a minimum  total  number  of  steps. 

3.1.3  Variables  of  Integration  for  This  Application 

The  array  of  variables  of  integration  X and  their  normaliza- 
tion values  X are  shown  in  Table  1.  The  numbers  identifying  the  variable 
locations  in  the  array  correspond  to  the  column  headed  by  NERR  in  the 
trajectory  time  history  output.  Thus  NERR  identifies  the  variable  having 
the  maximum  normalized  truncation  error. 


TABLE  1 VARIABLES  OF  INTEGRATION 


Location 

In  X Array 

Variable 

Definition 

Normalization 
Value  X 

1 

t 

time  (sec) 

1,000. 0 

2 

flight  path  angle  (rad) 

1.  0 

3 

V 

velocity  (ft/ sec) 

5,000. 0 

4 

Jk 

altitude  (ft) 

100,000.  0 

5 

R 

range  (ft) 

1,000,000. 0 

6 

W 

mass  (Ibm) 

1,000.  0 

7 

Ad 

ideal  velocity  (ft/sec) 

1,000, 0 

3.2  NEWTON-RAPHSON  PROCEDURE 

A Newton-Raphson  procedure  is  used  in  several  instances 
to  solve  a non-linear  equation  for  the  independent  variable.  Specifically, 
the  Newton-Raphson  method  is  used  to; 

(a)  Find  the  time  at  which  a phase  should  terminate  when 
a dependent  variable  stop  is  employed. 

( RUNGE  K subroutine) 

(b)  Find  the  angle  of  attack  required  to  satisfy  a given 
path  angle  rate.  ( XALPHl  subroutine) 
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(c)  Find  the  angle  of  attack  required  to  satisfy  a given  normal 
acceleration,  ( XALPH2  subroutine) 

(d)  Find  the  angle  of  attack  required  to  satisfy  a path  angle 
rate  which  varies  during  the  iteration, 

(DERIV  subroutine,  ICONT  = 2 option) 

Let  the  independent  variable  be  denoted  by  JC  , and  let 
the  non-linear  function  of  X be  , The  condition  to  be  satisfied  is 

^(JL)  s s 0 

where  is  the  desired  value  of  Aj  , The  derivative  of  is  found  by 

finite  differences  to  be 

fit)  = 

The  next  guess  for  X is 

“ Xi* 

, f 

The  process  is  repeated  until  X is  within  a given  tolerance  of  zero. 
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4.  SUBMODEL  STRUCTURE 


The  principal  subroutines,  their  primary  functions,  and 
the  flow  between  them  are  illustrated  in  Figure  3.  All  subroutines  are 
not  shown;  e.g.  , subroutines  XALPHl  and  XALPH2  used  for  angle  of  attack 
iterative  solutions  by  DERIV,  subroutines  used  in  support  of  PROPl,  the 
ramjet  performance  subroutine,  and  FORTRAN  supplied  subroutines  are 
omitted. 

All  subroutines  used  which  are  not  FORTRAN  supplied 
are  listed  and  briefly  described  in  Table  2. 
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TABLE  2 VEHPER  SUBROUTINES 


Trajectory  Related  Subroutines 

NAME 

DESCRIPTION 

AIRl 

Initializes  atmosphere  computations 

AIR 

Gets  atmospheric  properties  on  trajectory  (an  AIRl 
entry  point) 

DERIV 

Computes  derivatives  of  integrat  on  variables 

DSLINE 

Linear  interpolation  for  functions  of  two  independent 
variables 

ERROUT 

Not  Used 

MAINS 

Initializes  phase 

OUTPUT 

Prints  trajectory  time  history 

RUNGEK 

Integrates  simultaneous  first-order  differential 
equations 

STDATA 

Initializes  trajectory 

SEINE 

Linear  interpolation  for  functions  of  a single  indepen- 
dent variable 

TPROP 

TPROP  has  four  entry  points.  TPROP  entry  point 
initializes  linearized  ramjet  model. 

TPROPA 

Entry  point  tests  flight  conditions  for  region  bounds 

TPROPl 

Entry  point  gets  maximum  or  minimum  thrust  engine 
performance 

TPROP2 

Entry  point  gets  throttleable  engine  performance 

XALPHl 

Uses  Newton- Raphs on  method  to  get  angle  of  attack 

XALPH2 

Finds  angle  of  attack  and  thrust  required  to  satisfy 
normal  and  tangential  accelerations 

VEHPER 

Starts  trajectory  com''utation  by  calling  STDATA 

TABLE  2 (Continued) 


NAME 

DTRGET 

ISEN 

MACHNO 
PROP 1 1 
RGAMER 

TLUl 

TLU2 


Ramiet  Propulsion  Subroutines 
DESCRIPTION 

Determine  temperature  rise  by  table  look-up 

Determines  real  gas  total  temperature  and  pressure 
from  static  temperature  and  pressure 

Determine  Mach  No.  from  area  ratio  change 

Principal  ramjet  performance  subroutine 

Determines  ratio  of  specific  heats  (gamma)  and  gas 
constant  from  table  look-up 

Performs  linear  interpolation  for  functions  of  a single 
independent  variable 

Performs  linear  interpolation  for  functions  of  two 
independent  variables 


Turbojet  Propulsion  Subroutines 

NAME 

DESCRIPTION 

TJPER 

Computes  turbojet  off-design  performance  (thrust  and 
fuel  flow)  as  a function  of  flight  conditions  (angle-of- 
attach,  altitude,  and  Mach  No.  ) and  thrust  required. 

The  principal  function  of  TJPER  is  to  match  inlet  air- 
flow and  pressure  recovery  to  the  engine  requirement. 

GENENG 

A package  of  subroutines  which  for  given  flight  condi- 
tions and  pressure  recovery,  computes  thrust,  fuel 
flow,  and  corrected  airflow 

INLET 

A package  of  subroutines  which  for  given  flight  condi- 
tions and  corrected  airflow,  computes  inlet  pressure 

recovery 

FIGURE  3 SUBMODEL  FUN"IONAL  FLOW 


int«grate  phase  TPROP,  TPROPA OUTPUT 

Linearise  ramjet  data  Compute  and  print 

trajectory /output 

data 
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